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This  documentation  effort  was  supported  by  the  United  States  Army 
Missile  Command  under  Contract  No.  DAAH01 -73-C-l 150,  titled  "TACOC  II 
Model  Documentation.” 

The  documentation  is  a  nine-part,  three-volume  report.  Volume  I  is 
the  EXECUTIVE  SUMMARY,  Volumes  i IA,  B,  C,  and  0  are  the  PROGRAMMER/ ANALYST 
MANUALS  for  FRAG1 ,  FRAG 2 ,  FRAG3 ,  and  TERAIN,  PMAP  and  S0RTEV,  respectively. 
Volumes  IIIA,  B,  C,  and  0  are  the  USER/PLANNER  MANUALS  for  FRAG1 ,  FRAG2, 

FRAG 3 ,  and  TERAIN,  PMAP  and  S0RTEV,  respectively. 

As  an  aid  to  the  user,  a  somewhat  different  form  of  paragraph  and 
page  numbering  scheme  has  been  incorporated  in  this  documentation.  Each 
chapter  title  is  the  name  of  that  portion  of  TACOS  being  discussed.  There 
Is  one  exception  to  this  rule:  the  first  chapter  of  each  of  Volumes  II 
and  III  is  a  general  overview  of  the  TACOS  operation.  It  is  titled  "Chapter 
TACOS."  The  headings  within  each  chapter  are  of  a  modified  numerical 
scheme.  Except  in  a  few  places,  the  numbering  is  held  to  four  levels.  For 
each  volume,  the  level  designation  and  accompanying  heading  naw?»s  typography 
is: 


Illustrations  and  figures  are  numbered  consecutively  within  each  section. 
The  chapter  and  section  number  are  an  integral  part  of  the  numbering  scheme 
(i.e.,  Figures  FIA.l-1,  FIA.1-2,  etc.).  Page  numbering  uses  essentially 
the  same  scheme  used  for  heading  designations,  however,  only  three  levels 
are  used.  The  level  indicated  in  the  page  number  corresponds  to  the  major 
chapter  division  level  on  a  page.  An  example  would  be: 
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F3C. 1.2-1,  2,  3  ...  n. 

The  user  needs  simply  to  locate  the  section  of  interest  in  the  Table  of 
Contents;  he  can  then  turn  directly  to  the  appropriate  page.  A  variation 
on  this  methodology  is  used  in  the  documentation  for  individual  programs. 
The  program  name  being  documented  or  flowcharted  is  imbedded  in  the 
pagination  .;heme  at  the  third  level.  In  the  documentation  portion,  pages 
are  numbered  consecutively  with  letters  of  the  alphabet  i.e., 

FIA.4.MA1N-A,  B,  C  ... 

In  the  flowchart  portion,  pages  are  numbered  consecutively  with  numbers 
corresponding  to  the  flowchart  page  i.e., 

F1A. 4. MAIN-1 ,  2,  3  ... 

This  numbering  scheme  is  somewhat  nonstandard,  but  it  is  designed  to  afford 
the  reader  maximum  ease  of  use. 

Principal  contributers  to  this  work  include:  D.  E.  Edgemon,  J.  N. 
Gant,  D.  R.  Jackson,  J.  S.  Nowick? ,  J.  J.  Sikora,  L.  H.  Skinner,  and 
R.  J.  Upham.  Project  leadership  was  by  R.  L.  Katz  under  the  directorship 
of  0.  V.  Fedor off. 
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CHAPTER:  PMAP 

PMAP  PROGRAMMER/ANALYST  GUIDE 

P.l  MODULE  STRUCTURE 

The  structure  of  the  PMAP  module  Is  defined  by  specifying 
the  functional  flow,  a  detailed  narrative  description  of  each  submodel 
contained  in  PMAP,  a  description  of  inputs  required  and  outputs  generated 
by  PMAP  and  the  subroutine  call  sequence. 


P.1.1 


Functional  Flow 


Figure  P.1-1  illustrates  the  overall  functional  flow  of 
module  PMAP.  As  shown,  PMAP  reads  the  input  control  cards,  checks  the  site 
data  as  the  sites  are  input,  plots  a  site  map,  checks  the  path  data  as  the 
paths  are  input  and  plots  the  input  paths. 


P.1.2 


P.1.3 


Overlay  Structure 


PMAP  is  not  overlaid. 
Subroutine  Call  Sequence 


The  PMAP  subroutine  tree  is  diagrammed  in  Figure  P.l -2 .  The  sub¬ 
routine  call  sequence  is  indicated  in  Figure  P.I-3.  The  name  of  each  called  sub¬ 
routine  is  indicated  in  the  upper  left-hand  corner  of  each  box.  "Nesting"  of  sub¬ 
routine  calls  is  indicated  by  the  nesting  of  boxes.  This  figure  shows  all  calls 
which  may  be  made  in  the  course  of  a  run.  This  does  not  include  calls  to  ERR  which 
is  called  only  when  an  error  is  detected.  Also,  the  diagram  does  not  include  caJIs 
to  any  system  supplied  utility  routines. 
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? MAP  Detailed  Narrative  Description 


The  objective  of  this  detailed  narrative  description  is  to 
provide  a  brief,  narrative,  non-mathematlcal ,  synopsis  of  each  submodel 
in  the  subject  module.  Additionally  a  secondary  objective  is  to  provide 
some  correlation  of  submodels  wl.'.xin  this  nodule  and  some  other  modules 

where  close  correlation  exists.  A  staienwni  -liy-stnlrmrni  Ipt  Ion  »l 

\ 

the  code  contained  iri  each  subroutine  may  be^Found  In  the  subroutine  docu¬ 
mentation  and  associated  flew  charts  in  Paragraph  PA. 
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Figure  P.l-1.  PMAP  Functional  Flow 
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BRADDOCK,  DUNN  AND  McDONALD. INC. 


C1TN 

(ENTRY  TO  C00RT)  SEE  BLOCK  2 
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0NEPT  SEE  BLOCK  I 


0NEPT  SEE  BLOCK  I 


0NEPT  SEE  BLOCK  i 


ZN0PLT 
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S0RTR  SEE  BLOCK  4 


S0RTR  SEE  BLOCK  4 


MAP  SEE  BLUCK  5 


Figure  P.I-3.  PMAP  Subroutine  Call  Sequence  (Continued) 
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BLOCK  I 


Figure  P.I-3.  PMAP  Subroutine  Cal!  Sequence  (Continued) 
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figure  P.I-3.  PMAP  Subroutine  Call  Sequence  (Continued) 
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P. 1.4.1  Sequence  of  Operation 

As  shewn  by  Figure  P.l-I  (PMAP  Functional  Flow),  the  sequence 
of  operation  is  basically  to  read  the  input  systems  and  site  cards,  perform 
the  specified  site  data  checks  and  site  plotting,  read  the  input  path  data, 
perform  the  specified  path  checks  and  piofc  the  paths,  and  determine  site/ 
oath  ei'yageabi  I  i  ty. 

P.1.4,2  Terrain  Data  Base 

PMAP  utilizes  the  identical  real-world  terrain  data  base  as 
does  FRAG1A  and  FRAG18.  This  data  base  is  used  to  check  for  terrain  avail¬ 
ability  at  the  site  locations  and  path  points  and  to  check  the  input  site 
altitudes . 

P.1.4. 3  Site  Data  Checking 

As  the  SITE  cards  are  input  to  PMAP,  they  are  checked  for 
redundancies,  such  as  duplicate  site  cards,  sites  with  duplicate  names 
and  sites  with  identical  locations.  If  requested,  each  site  is  also  checked 
for  terrain  data  related  errors.  Anomalies  are  entered  in  the  site 
description  table. 

P.1.4. 4  Path  Oata  Checking 

As  the  path  data  are  read,  each  path  is  checked  to  determine 
if  it  contains  an  unusually  long  path  !eg  and  each  turning  point  is  checked 
for  an  excessive  turn  angle.  Anomalies  are  entered  in  the  path  description 
table. 

P.1.4.5  Targeted  Path  Identification 

PMAP  identifies  paths  which  are  targeted  against  specific 
sites  by  matching  analyst  specified  path  point  identifiers  with  site  identi¬ 
fiers.  targeted  path  points  are  identified  in  the  path  description  table. 


P.l-U 
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P. 1.4.5  Path  Enqageabi ! i ty  Identification 

For  each  site,  PMAP  determines  the  paths  which  pass  within 
the  site  range  of  interest.  These  paths  are  considered  engageable  by  this 
site  ano  are  entered  for  that  site  in  the  site  engageabi I i ty  table. 

P. 1 .4.6  Hap  Piott ing 

A  PMAP  plot  consists  of  printing  the  site  identifiers  at 
the  site  locations  and  the  path  point  sequence  numbers  at  the  point  locations 
on  a  UTm  grid.  One  plot  is  printed  for  aP  trie  site  locations  and  one  plot 
for  each  of  the  input  paths. 

P.2  INPUT  RECORD  FORMATS 

P.2.1  Digitized  Terrain  Data  Tape 

The  digitized  terrain  data  tape  utilized  f  /  PMAP  has  the 
identical  format  as  the  terrain  data  tape  utilized  by  FRaGIA.  Fo>-  the 
detai  led  tape  format  see  F!A.<!,  Volume  II  A  FRAG1  Prog  rammer /Ana  I  yst  Manual. 

P . 2 . 2  Control  Cards 

In  order  to  perform  the  preliminary  data  checks  and  plot 
the  site  and  path  .saps,  PMAP  expects  the  following  parameter  types: 

•  Terrain  file  specification 
«  Control  logicals 

•  Systems  data 

•  Site  specifications 

•  Path  specifications 

A  more  detailed  explanation  of  the  card  type  inputs  is  found 
in  Volume  111  D  TERRAIN,  PMAP,  S0RTEV  User/Pla.iner  Marual,  Paragraph  P.2. 

P.3  MATHEMATICAL/-. 03 !CAl  STRUCTURE 

PMAP  utilizes  the  following  five  submodels:  Site  Data  Check¬ 
ing,  Path  Data  Checking,  Targeted  Path  Identification,  Path  Engageabi 1 i ty 
Identification,  and  Map  Plotting. 


P.2-1 
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P.3.1  Site  Data  Checking 

P.3.1.1  Purpose; 

Tnis  submodel  Identifies  probable  errors  in  site  specifica¬ 
tion. 

P.3. 1*2  Aigorithm/Function  Definition; 

Checks  are  perfoimed  on  a  card  by  card  basis. 

P.3.1.3  Rationale  for  Derivation: 

As  each  SITE  card  is  read,  the  terrain  data  directory  Is 
accessed  to  determine  if  the  terrain  file  contains  the  grid  square  con¬ 
taining  the  site  location.  If  terrain  data  are  available  for  this  site 
location,  the  terrain  c'n„.  are  used  to  interpolate  for  the  terrain  elevation 
at  the  site  location.  If  terrain  is  not  available  at  the  site  location 
or  the  site  altitude  specified  is  less  than  the  terrain  elevation,  an 
informative  message  is  printed  in  the  site  description  table. 

The  site  elevation  is  determined  by  bilinear  interpolation 
of  the  four  closest  terrain  data  sample  points.  The  procedure  is  as  fol¬ 
lows: 


4 

(x, ,y2) 

a  - -  -  ■ 

(x2>2) 

XREM*RATE 

* 

YREM*RAT£ 
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-  -- 

>  — - 

>~  ^ 

>< 

” — - * - 

'  iX2,Y1 ) 
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SITE  LOCATION 
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Where:  RATE  is  the  lateral  data  rate 
[k]  is  ’’the  integer  part  cf  K" 
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A  card  by  card  check  Is  than  performed  to  determine  If  a 
site  card  has  been  previously  processed  whose  site  Identifier,  location, 
or  both  Identifier  and  location  are  identical  with  the  card  under  cons  I dera 
tlon.  If  any  of  these  duplications  are  detected,  an  "error''  message  is 
printed  in  the  site  description  table. 

P.3.1 Jf 
P.3. 1.5 
P.3.1.6 
P.3.2 
P.3.2.1 

tlon. 

P.3.2.2 

P.3.2.3  Rationale  for  Derivation 

As  each  'P0INT1  card  is  read,  the  two  letter  grid  square 
designator  for  the  point  location  is  matched  against  the  valid  designators 
from  the  terrain  file  directory  to  determine  If  this  point  lies  on  avail¬ 
able  terrain.  Upon  entry  of  the  second  and  ail  succeeding  cards,  the 
leg  length  of  the  preceding  leg  is  calculated  and: this  length  check¬ 
ed  to  determine  if  it  is  longer  than  the  maximum  specified.  Upon  entry  of 
the  third  and  all  succeeding  points,  the  path  turn  angle  Is  tested  to 
dtMrrminr  if  (he  maximum  1  urn  angle  fh**  pi'rviouvly  J tiptit  paint  h.v». 
exceeded  the  maximum  expected  turn  angle.  The  test  for  maximum  turn  angle 
is  developed  as  follows: 


P.3-3 
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TRACK  1 
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TRACK  1 
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TRACK  2 
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/-BLOCK  POINTER  FOR  TRACK  1 
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Figure  P.3-1  Path  Point  Storage 
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P.3.2.4  Parameter  Definitions:  None 

P.3.2.5  Assumptions  and  Limitations:  None 

P.3.2.6  Data  Edit  Regui rements  and  Manual  Procedures:.  None 

P.3.3  Targeted  Path  Identification 

P.3.3.1  Purpose 

This  identifies  explicitly  the  path  points  which  are  target¬ 
ed  for  a  specified  site. 

P.3.3.2  Algorithm/Function  Definition 

A  path  is  considered  to  be  targeted  for  a  site  if  an  input 
path  point  location  is  within  25  meters  ground  range  cf  that  site  location. 

P.3.3.3  Rationale  for  Derivations 

As  each  path  point  is  read,  the  point  identifiers  are  compared 
with  each  of  the  site  identifiers  to  determine  if  this  point  is  a  candidate  for 
targeting  against  any  of  the  stored  sites.  If  a  match  is  determined,  the 
range  from  that  site  to  the  path  point  is  calculated.  This  range  is  then 
tested  and  if  it  is  less  than  25  meters,  this  path  is  targeted  against  this 
site. 

P.3.3.4  Parameter  Definitions:  None 

P.3.3-5  Assumptions  and  Rationale:  None 

P.3.3.6  Data  Edit  Requirements  and  Manual  Procedures: 

If  a  path  point  is  to  be  considered  for  targeting,  its 
Identifier  must  match  the  site  identifier  under  consideration. 

P - 3 • 4  Path  Engageabi I i  iy  Idem  i  ficat  inn 

P.3.4.1  Purpose 

This  submodel  determines  which  paths  intersect  each  site 
engagement  volume. 

zo< 
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P.3.4.2  Algorithm/Function  Definitions 

Whenever  a  vehicle  passes  within  a  site's  range  of  interest, 
it  is  considered  engageable. 

P.3.4.3  Rationale  for  Derivation 

For  each  site/path  combination,  the  vehicle  path  is  checked, 
path  leg  by  path  leg,  to  determine  if  this  path  leg  intersects  this  site's 
engagement  volume.  Each  path  leg  is  tested  utilizing  the  rationa’e  developed 
in  Figure  P.3“2.  First,  the  square  of  the  distance  between  the  path  leg 
end  points  (8  ),  the  square  of  the  range  from  the  site  to  path  point  i  (C  ) 
and  the  square  of  the  range  from  the  site  to  the  path  point  i+l  (A  )  are 
calculated  by  summing  the  squares  of  the  differential  easting,  northing  and 
altitudes.  The  cosine  of  the  angle  y  is  then  determined  by  the  law  of  cosines. 
The  side  8'  of  triangle  A,  B',  R'  is  then  determined  as  shown.  If  it  is  then 
determined  that  R'  is  less  than  the  range  of  interest,  the  path  intersects 
this  site's  engagement  volume. 

P.3.4. 4  Parameter  Def ini t ions: 

P.3. A. 4.1  Range  of  interest 

That  range  within  which  it  is  possible  for  an  air  defense 
site  to  detect  an  airborne  vehicle. 

P.3.4.4.1.1  Source:  System  Description 

P.3.4.4.1.2  Representative  Values: 

35000,78000  meters. 


P.3.4.4.1.3 

P.3.4.5 
P.3. 4.6 


Ranges : 

1000,000  meters 

Assumptions  and  Rationale:  None 

Data  Edit  Requirements  and  Manual  Procedures:  None 
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P.3.5  Hap  Plotting 

P.3.5.1  Purpose: 

Provide  graphical  displays  of  site  locations  and  vehicle 

paths. 

P.3.5.2  Algorithm/Function  Definitions: 

PHAP  plots  consists  of  printing  the  site  identifiers  at  the 
site  locations  and  the  path  point  sequence  numbers  at  the  point  locations 
on  a  UTM  grid. 

P.3.5.3  Rationale  for  Derivation: 

Once  the  points  to  be  plotted,  either  site  locations  or  path 
points,  have  been  determined,  the  plot  algorithm  begins  by  first  sorting  the 
points  and  restoring  them  in  the  order  of  descending  northing.  Next,  the 
northern-most,  southern-most  and  western-most  points  are  determined.  The 
map  plotting  algorithm  then  uses  the  northern-most  point  to  set  the  map 
northern  boundary,  the  western-most  point  to  set  the  weastern  boundary  and 
the  southern  most-point  to  set  the  southern  boundary.  The  width  of  the 
strips  to  be  plotted  is  then  determined  from  the  plot  scale  and  page  width. 

The  easting  of  the  eastern  boundary  of  the  first  strip  to  be  plotted  is  then 
dete. mined  by  subtracting  the  strip  width  from  the  easting  of  the  western 
boundary.  This  north-south  strip  is  then  printed,  line  by  line,  with  the 
UTM  grid  and  point  identifiers  printed  as  required.  If  all  the  points 
under  consideration  did  not  lie  on  this  strip,  the  algorithm  mov.is  one 
strip  east  and  prints  the  next  north-south  strip.  This  process  continues 
until  all  the  points  under  consideration  have  been  plotted. 

p.3.5.4  Parameter  Definitions :  None 

p.3.5,5  Assumptions  and  Rationale:  None 

p.3.5  $  Data  Edit  Requirements  and  Manual  Procedures:  None 
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P.4  PROGRAM  LOGIC 

P . 4 . 1  Chapter  Arrangement 

The  following  chapter  presents  detailed  program  Information 
about  the  routines  which  comprise  PHAP. 

The  routines  are  presented  in  alphabetical  order  arranged 
in  two  parts.  For  each  routine  the  first  part  consists  of  a  13  point 
documentation  and  the  second  part  consists  of  an  annotated  logic  f  low- 
chart - 

P.4.2  Flowchart  Conventions 

Off-page  connectors  are  of  two  types:  (1)  to  following 
page  and,  (2)  to  another  page.  A  Type  1  connector  occurs  only  at  the  end 
of  a  page  and  contains  the  number  of  the  page  to  which  it  is  connected. 

The  same  number  appears  in  a  connector  box  at  the  top  of  its  associated 
page.  The  Type  2  connector  consists  of  two  boxes.  The  upper  box  con¬ 
tains  the  page  number  to  which  the  connector  points.  The  lower  box 
contains  the  alphanumeric  designator,  on  that  page,  to  which  the  con¬ 
nector  points,  rcr  example,  on  Page  2  of  subroutine  T0DD  in  S0R1CV, 
appears  an  off-page  connector  labeled  ^ .  This  means  that  this  con¬ 
nector  points  to  a  connector  labeled  A1  which  is  on  Page  3.  Additionally, 
the  1  Indicates  it  is  the  first  connector  referenced.  Halfway  dowr.  Page 
3  is  another  connector  labeled  A2.  These  connectors  are  labeled  A! 
through  A9,  Bl  through  B9,  etc.  Additionally,  each  Type  2  on-page  connec¬ 
tor  is  accompanied  by  a  label  of  the  page  (or  pages)  number  that  references 
this  connector.  Connections  between  the  first  and  last  statements  of  a 
DO  loop  are  indicated  by  a  dashed  line.  If  these  statements  are  on 
separate  pages,  the  dashed  line  from  the  first  statement  terminates  with 
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a  pointer  to  the  appropriated  page.  The  dashed  line  for  the  last  state¬ 
ment  on  the  other  page  terminates  with  a  pointer  back  to  the  originating 
page. 

Each  FORTRAN  statement  number  appears  to  the  upper  left 
of  its  appropriate  logic  flow  box. 

The  subroutines  comprising  PMAP  are: 


ABSC0R 

MAXN 

S0RTR 

CLN 

MVC 

STPNT 

C00RT 

MYTIME 

THYME 

CVNBM 

0NEPT 

TT0RNG 

EBCO 

PACK 

TTRK 

ELEV 

PMAP 

ZCKPTJ 

POO 

PAW 

JTRK 

SEARCH 

LABEL 

SET8PT 

MAIN 

MAP 
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1 .  NAME :  ABSC0R 

2.  TYPE  Or  PROGRAM:  SUBROUTINE,  LOGICAL 

3.  LANGUAGE:  FORTRAN  EXTENDED 

k.  PURPOSE:  Fill  the  print  buffer  with  the  abscissa  coordinate 
characters  to  be  printed. 

5.  ASSUMPTIONS  AND  LIMITATIONS:  The  number  of  labels  must  fit  on 
the  page  at  the  selected  scale  factor. 

6.  ERROR  RETURNS:  A  call  ERR  and  a  RETURN  A1  will  be  executed  if 
the  following  errors  are  detected: 

A.  The  abscissa  coordinate  values  exceed  132  spaces  at  the 
specified  spacing. 

B.  An  illegal  grid  letter  was  returned  by  CNL. 

7-  LINKAGE  STATEMENT  AND  DESCRIPTION  OF  ARGUMENTS: 

CALL  ABSC0R{ I  START,  ISPC,  N) ,  RETURNS  (Al) 

ISTART  -  Integer,  the  starting  position  for  the  abscissa  coordinates 
in  the  print  array. 

ISPC  -  The  number  of  blanks  between  coordinate  labels. 

N  -  Number  of  labels  to  be  printed  across  the  page. 

Al  -  Exit  taken  if  one  of  the  conditions  listed  under  6  (above) 
occurs . 

8.  PROGRAM  CALLING  THIS  PROGRAM:  MAP 

9.  COMMON  VARIABLES  USED: 

A.  ZZC0RD  -  DELTA,  ESTEDG,  I0RBUF,  WSTEDG 

B.  ZZCVi  -  X0RG 
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10.  COMMON  VARIABLES  TO  BE  SET: 

A.  ZZC0RD  -  DELTA,  ESTEOG  ,  WSTEDG , 

B.  ZZCVI  -  X0RG 

11.  COMMON  VARIABLES  CHANGED: 

A.  ZZC0RD  -  I0RBUF 

12.  PROGRAMS  USED  AND  DESCRIPTION  OF  LINKAGES: 

A.  SUBROUTINES 

(1)  CALL  ERR  (NAME,  LAST,  NSTAT)  -  Prints  an  error  message 
indicating  name  of  subroutine  and  statement  number 
where  error  occurred. 

NAME  -  Name  of  subroutine  where  error  was  detected. 

LAST  -  End  address  of  the  storage  area  to  be  dumped  from 
subroutine  NAME. 

NSTAT  -  Statement  number  in  subroutine  NAME  where  error 
occurred . 

(2)  CALL  LABEL  (M,  X,  X0RG,  ID),  RETURNS  (Al) 

M  -  Character,  contains  the  one  letter,  three  digit 
grid  ’abel  for  the  line  X. 

X  -  Relative  easting  (ID*E)  or  northing  (|D«N)  for 
which  the  UTM  grid  letter  is  desired. 

X0RG  -  Displacement  (meters)  of  the  origin  from  the  UTM 
reference  point  (AA) . 

ID  -  Character,  E  indicates  easting  of  strip  is  desired 
N  indicates  northing  of  line  is  desired. 

Ai  -  Control  is  transferred  to  this  statement  when  an 
illegal  grid  letter  is  detected  by  CLN. 
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(3)  CALL  MVC  (ICHT,  1NPTR,  ZINADD,  I0PTR,  ZI0ADD)  -  Hove 

characters  from  the  word  ZINADD  into  the  word  ZI0AOD. 

ICHT  -  Integer,  number  of  characters  to  be  moved  from 
ZINADD  to  ZI0ADD. 

INPTR  -  Integer,  index  of  the  first  character  to  be 
moved  from  ZINADD. 

ZINADD  -  Character,  2  word  array  (20  character)  from  which 
characters  will  be  moved. 

I0PTR  -  Integer,  index  of  the  location  in  ZI0ADD  where 
the  first  character  will  be  moved. 

ZI0ADD  -  Character,  2  word  array  (20  character)  f^om  which 
the  characters  will  be  moved. 

13.  NOTES  ON  METHODOLOGY:  Upon  entry  to  ABSC0R,  0PTR(l)  Is  set  to 
the  easting  of  the  western  edge  (WSTEDG)  of  the  map  and  M  is  set 
to  the  number  of  labels  to  be  printed  plus  one.  A  DO  loop  is  then 
entered  to  generate  the  coordinate  abscissa  labels  (1  letter, 

3  digits).  LABEL  is  used  to  generate  the  labels  (I0RPT)  at  DELTA 
meters  (intervals),  starting  at  the  western  edge.  This  loop 
generates  labels  until  the  eastern  edge  (ESVEDG)  of  map  is  passed 
or  the  maximum  number  of  labels  (N)  have  been  found.  If  N  labels 
are  generated  before  the  eastern  edge  is  passed,  J,  the  number  of 
labels  actually  generated,  is  set  to  N  and  control  transfers  to 
statement  1700,  If  the  eastern  edge  is  passed,  control 
transfers  to  statement  1100  where  the  last  label  is  set  to  a 
field  of  blanks.  At  statement  1700  the  starting  position  (ISTART) 
in  the  print  array  and  the  last  position  in  the  print  array  are 
tested.  If  the  start  position  is  less  than  one  or  the  last  position 
is  greater  than  132,  control  is  transferred  to  statement  9000  for 
error  processing. 


*  * 
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M  is  then  set  to  J  minus  one  and  a  DO  loop  executed  which 
places  blank  fill  in  all  the  entries  of  the  array  I0RBUF.  A  major 
DO  loop  is  then  entered  which  moves  the  k  characters  (one  tetter, 

3  digit)  representing  the  abscissa  coordinate  label  into  the  buffer 
I0RBUF.  The  counters  M0RE  and  LC  keep  track  of  the  location  in 
I0RBJF  where  the  characters  are  to  be  stored.  An  inner  DO  loop 
moves  the  characters  for  each  abscissa  into  the  buffer  I0RBUF. 

When  all  four  characters  for  each  abscissa  coordinate  label  have 
been  stored  in  I0RBUF,  control  transfers  to  the  calling  program. 

Upon  detection  of  an  error  in  A8SC0R,  control  is  transferred 
to  the  appropriate  error  statement  where  a  call  is  made  to  ERR, 
which  prints  an  error  message.  A  RFTURN  Al  is  executed  and  control 
transferred  to  the  calling  program. 
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1.  NAME:  .CLN  (other  entry  points:  iL) 

2.  TYPE  OF  PROGRAH:  SUBROUTINE,  LOGICAL 

3.  LANGUAGE:  FORTRAN  EXTENDED 

k.  PURPOSE:  Calculate  the  number  of  100  km  squares  between  a  reference 
point  and  given  grid  square  designator.  Given  the  number  of  squares 
between  the  reference  and  a  designator,  find  the  appropriate  grid 
square  designator. 

5.  ASSUMPTIONS  AND  LIMITATIONS:  None 

6.  ERROR  RETURNS:  A  RETURN  Al  will  be  executed  if  an  illegal  grid 
square  designator  is  found. 

7.  LINKAGE  STATEMENT  AND  DESCRIPTION  OF  ARGUMENTS: 

CALL  CLN  (LL,  N,  D) ,  RETURNS  (Al) 

entry  point  CALL  CNL  (LL,  N,  D) ,  RETURNS  (Al) 

LL  -  integer,  two  letter  grid  square  designator. 

N  -  Integer,  number  of  100  x  100  km  grid  squares  between  input  point 
and  the  reference  point  (location  AA). 

D  -  Indicates  whether  the  squares  counted  are  easting  (E)  or  northing 
(N). 

Al  -  Exit  is  taken  when  an  illegal  grid  square  designator  is  detected. 

8.  PROGRAMS  CALLING  THIS  PROGRAM: 

Entry  point  CLN:  C02RT 
Entry  point  CNL:  C00R7,  LABEL 

9.  COMMON  VARIABLES  USED:  None 

10,  COMMON  VARIABLES  SET:  None 

11.  C0HHQN  VARIABLES  CHANGED:  None 
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12.  PROGRAMS  USED  AND  DESCRIPTIONS  OF  LINKAGES: 

A.  SUBROUTINES 

(1)  CALL  MVC  (ICNT,  INPTR*  Z INADD,  10PTR,  ZI0ADD)  -  Move 

characters  from  the  word  ZINADD  into  the  word  ZI0ADD. 

ICNT  -  Integer,  number  of  characters  to  be  moved  from 
ZINADD  to  ZI0ADD. 

INPTR  -  Integer,  index  of  the  first  character  to  be 
moved  from  ZINADD. 

ZINADD  -  Character,  2  word  array  (20  character)  from  which 
characters  will  be  moved. 

I0PTR  -  Integer,  index  of  the  location  in  ZI0ADD  where 
the  first  character  will  be  moved. 

ZI0ADD  -  Character,  2  word  array  (20  character)  from  which 
the  characters  will  be  moved 

13.  NOTES  ON  METHODOLOGY:  CLN  provides  the  logic  necessary  to  count  the 
number  of  100  km  grid  squares  between  a  reference  location  AA  and  a 
specified  grid  square  corner.  By  ENTRY  CNL  the  one  letter  grid 
square  designator,  can  be  determined  from  a  count  of  the  number  of 
squares  between  the  reference  point  and  the  grid  square  corner. 

Upon  entry, GLN  sets  L  equai  to  the  one  letter  grid  square  desig¬ 
nator  of  the  reference  point  (A).  Next  MVC  is  called.  MVC  sets  L 
equal  to  the  first  letter  of  the  grid  square  designator  LL.  I,  the 
number  of  legal  grid  square  letters  is  next  determined.  If  D  is  equal 
I  HE  *  an  easting  count  is  desired  and  I  is  set-to  the; number  of  legal 
easting  points  [2k),  If  D  is  equal  INN,  a  northing  count  is  desired 
and  I  is  set  f.to  the  number  of  legal  northing  points  (20).  A  DO  loop 
tests  for  a  legal  grid  square  designator  from  the  DATA  table  LT  until 
a  grid  square  designator  is  found  which  equals  the  grid  square  desig¬ 
nator  L.  When  this  equalicy  if  found,  control  is  transferred  to 
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statement  1115  where  the  number  of  squares  (N)  is  determined  to  be 
the  loop  counter  minus  one;  A  normal  RETURN  is  executed  and  control 
transferred  to  the  calling  subroutine.  I f  no  match  is  found  for  the 
one  letter  grid  square  designator,  control  Is  transferred  to  the 
caj ling  subroutine  via  error  return  RETURN  Al. 

Upon  entry  at  ENTRY  CNL  a  test  is  made  on  the  square  counter  N. 

If  the  square  counter  (N)  is  less  than  zero  a  RETURN  A!  is  taken. 

If  N  is  not  less  than  zero,  I,  the  number  of  legal  grid  square 
designators,  is  set  to  23  if  easting  is  being  determined,  and 
set  to  19  if  northing  is  being  determined.  Next  the  dispatching 
number,  Nl,  is  determined.  Then  the  letter  for  the  desired  grid 
square  designator  LL  is  determined  from  the  dispatch  table  LT.  A 
normal  exit  is  taken  and  control  transferred  to  the  calling  subroutine. 
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1.  NAME:  C00RT  (other  entry  points:  CVLN,  CVNB,  CITA,  CITN,  SET0RG) 

2.  TYPE  OF  PROGRAM:  SUBROUTINE,  MATHEMATICAL 

3.  LANGUAGE:  F0RTRAN  EXTENDED 

k.  PURPOSE:  To  calculate  the  data  required  for  definition  of  a 
coordinate  system  origin  and  facilitate  the  conversion  from 
relative  easting  and  northing  to  UTM  and  UTM  to  relative  easting 
and  northing  for  this  system. 

5.  ASSUMPTIONS  AND  LIMITATIONS:  The  data  required  for  conversion  of 
UTM  to  relative  easting  and  northing  coordinate  must  be  entered 
via  the  argument  list.  AH  other  conversions  must  have  the  appro¬ 
priate  data  in  common. 

6.  ERROR  RETURNS:  A  RETURN  A I  will  be  executed  if  either  of  the 
following  errors  are  detected: 

A.  An  illegal  2  letter  UTM  grid  square  designator  has  been  found. 

B.  Point  being  converted  from  UTM  to  relative  easting  and  northing  lies 
west  or  south  of  the  coordinate  system  origin. 

7.  LINKAGE  STATEMENT  AND  DESCRIPTION  OF  ARGUMENTS: 
entry  point  CALL  CVNB  (X,  Y,  L) 

entry  point  CALL  CITA  (X,Y,  L),  RETURNS  (Al) 
entry  point  CALL  C1TN  (X,  Y,  L)t  RETURNS  (Al) 
entry  point  CALL  SETUPS  (X,  Y,  L) ,  RETURNS  (Al) 

X  -  Relative  easting  of  the  point  under  consideration. 

Y  -  Relative  northing  of  the  point  under  consideration 
L  -  Two  letter  UTM  grid  square  designator. 

Al  -  Control  is  transferred  to  this  statement  number  if  an  error  is 
detected  in  C00RT. 
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8.  PROGRAMS'?  CALL  I NG  THI  S  PROGRAM: 

Entry  point  CVNB:  MAIN 
Entry  point  CITA;  MAIN,  0NEPT 
Entry  point  C1TN:  MAIN,  0NEPT 
Entry  point  SET0RG:  MAIN 

9.  COMMON  VARIABLES  USEO: 

A.  ZZCVI  •*  1 0RG ,  J0RG ,  LT0RG,  X0RG,  Y0RG 

B.  ZZCV2  -  LTRS ,  NBRS,  NBRBCD 

10.  COMMON  VARIABLES  TO  BE  SET: 

A.  ZZCVI  -  I0RG,  J0RG,  LT0RG,  X0RG,  Y0RG 

B.  ZZCV2  -  LTRS,  NBRS 

11.  COMMON  VARIABLES  CHANGED: 

A.  ZZCVI  -  I0RG,  J0RG,  LT0RG,  X0RG,  Y0RG 

B.  ZZCV2  -  LTRS,  NBRS,  NBRBCD 

12.  PROGRAMS  USED  AND  DESCRIPTION  OF  LINKAGES: 

A.  SUBROUTINES 

(I)  CALL  CLN  (LL,  N,  D) ,  RETURNS  (AI)  -  Calculates  the  number 
of  100  km  grid  squares  between  the  reference  point  AA  and 
and  specified  grid  square  designator. 

LL  -  Two  letter  UTM  grid  square  designator. 

N  -  Integer,  number  of  100  km  grid  squares  between  the 
specified  point  and  the  reference  point  \AA). 

D  -  Indicates  whether  the  squares  to  be  counted  are  east¬ 
ing  (E)  or  northing  (N). 
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A1  -  Exit  Is  taken  when  an  illegal  grid  square  designator  Is 
detected  in  CLN. 

(2)  CALL  CNL  (LL,  N,  0) ,  RETURNS  (Al)  (entry  point  to  CLN)  - 
Determines  the  corresponding  UTM  grid  square  designator  a 
specified  number  of  100  km  grid  squares  from  the  reference 
point  AA. 

LL  -  One  letter  grid  square  designator  corresponding  to 
either  the  northing  or  easting  direction  from  the 
reference  location  AA. 

N  -  Integer,  number  of  100  km  grid  squares  between  the 
•specified  point  and  the  reference  point  (AA) . 

D  -  Indicates  whether  the  squares  between  the  reference 
point  and  calculated  grid  square  designator  correspond 
to  an  easting  or  northing  direction. 

Al  -  Exit  is  taken  when  an  illegal  grid  square  designator  is 
detected  in  CNL. 

(3)  CALL  EBCD  (NBRS ,  IDUH,  NBRBUF)  -  Converts  the  8  digit  UTM 
grid  coordinate  into  an  alphanumeric  code  suitable  for 
printing. 

NBRS  -  8  digit  UTM  grid  coordinate. 

IDUH  -  Dummy  variables  not  used. 

NBRBUF  -  Array  of  dimension  3,  NBRBUF  (2)  containing  the 

coordinate  corresponding  to  easting  and  NBRBUF  (2) 
containing  the  coordinate  corresponding  to  northing. 

(k)  CALL  MVC  (ICNT,  INPTR,  Z INADD,  I0PTR,  ZI0ADD)  -  Move 
characters  from  the  word  ZINADD  into  the  word  ZI0ADD. 

ICNT  -  integer,  number  of  characters  to  be  moved  from 
ZINADD  to  ZI0ADD. 
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INPTR  -  Integer,  index  of  the  first  character  to  be 
moved  from  Z INADD. 

ZINADD  -  Character,  2  word  array  (20  character)  from  which 
characters  will  be  moved. 

I0PTR  -  Integer,  index  of  the  location  in  ZI0ADD  where 
the  first  character  will  be  moved. 

ZI0ADD  -  Character,  2  word  array  (20  character)  from  which 
the  characters  will  be  moved. 

13*  NOTES  ON  METHODOLOGY:  Entry  to  C00RT  via  entry  point  SET0RG  calcu¬ 
lates  the  displacement  of  the  origin  (east  and  north)  from  the  UTM 
reference  point  AA.  Upon  entry  to  SET0RG,  the  two  letter  UTM  grid 
square  designator  of  the  origin  (L)  is  stored  in  the  origin  data 
table  ZZCV1  and  the  UTM  data  table  ZZCV2,  and  logical  ZH0ME  is  set 
to  TRUE  to  insure  origin  setup.  At  statement  1151,  CLN  is  called 
to  find  how  many  100  km  grid  squares  east  of  the  reference  point 
AA  the  coordinate  system  origin  is  located.  The  displacement 
of  the  origin  in  meters  east  of  AA  is  then  calculated  and  stored  in 
integer  format  and  real  format  in  local  variables  II  and  XI,  respec¬ 
tively.  MVC  and  CLN  are  used  to  find  how  many  squares  north  of 
the  reference  point  the  coordinate  system  origin  is  located.  MVC  moves 
the  second  letter  of  the  UTM  grid  square  designator  into  the  first 
position.  CLN  finds  the  number  of  squares  north  of  the  reference 
point  the  origin  is  located.  The  displacement  of  the  origin 
in  meters  north  of  the  reference  point  is  then  calculated  and 
stored  in  integer  format,  and  real  format  in  local  variables  Jl 
and  YI,  respectively,  if  ZH0ME  is  FALSE,  the  displacement 
of  the  origin  is  stored  in  the  origin  data  table  and  control  is 
returned  to  the  calling  program. 
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If  entry  to  C00RT  Is  via  entr.y  CVLN  the  displacement  of  the  2 
letter  grid  square  designator  LTRS  from  the  reference  point  AA  Is 
determined.  Entry  CVLN  sets  logical  ZH0ME  to  TRUE  and  uses  the  same 
code  as  entry  SET0RG  to  determine  the  east  and  north  displacements 
from  AA.  These  displacements  (east  and  north)  are  stored  In  local 
variables  in  integer  format'  and  real  format  and  control  Is  returned 
to  the  calling  program. 

Entry  via  C1TN  converts  UTH  coordinates  to  relative  easting  and 
northing  for  the  coordinate  system  establ ished  via  entry  SET0RG. 

CITN  first  calculates  the  easting  and  northing  due  to  the  8  digit  UTH 
grid  coordinates  (NBRS).  The  call  to  CLN  at  statement  1171 
determines  how  many  grid  squares  east  of  the  reference  point  AA  the 
grid  square  LTRS  is  located.  A  call  to  HVC  and  CLN  finds  how  many 
grid  squares  north  of  the  reference  point  AA  the  grid  square  LTRS  is 
located.  The  relative  easting  (X)  and  northing  (Y)  of  the  point  of 
interest  is  then  determined  by  addition  of  these  components.  If  the 
point  lies  west  of  the  coordinate  system  origin  the  error  exit 
RETURN  Al  is  taken.  Otherwise  a  normal  return  is  taken  to  the  calling 
subroutine. 

Entry  CITA  converts  the  relative  easting  and  northing  of  a 
prespecified  point  to  UTH  coordinates.  Upon  entry  via  CITA, LTRS 
is  set  to  a  field  of  blanks,  I  is  set  to  the  number  of  10  meter 
intervals  between  the  specified  point  and  the  reference  point 
measured  along  the  east  axis,  and  k  is  found  to  be  the  number  of 
ICO  km  squares  between  the  reference  point  and  the  specified  point 
measured  along  the  east  axis.  At  statement  1 1 1 1,  CNL  is  called  to 
find  the  easting  letter  for  the  grid  square  designator.  Then  J 
is  set  to  the  number  of  iO  meter  intervals  between  the  specified 
point  and  the  reference  point  measured  along  the  north  axis  and 
K  is  found  to  be  the  number  of  100  km  squares  between  the  reference 
point  and  the  speci fied  point  measured  along  the  north  axis.  CNL  is 
called  at  statement  1112  to  find  the  northing  letter  corresponding 
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to  this  square.  HVC  Is  then  called  to  mask  together  the 
easting  letter  (LTRS)  and  the.  northing  letters  (H)  to  form  the  UTM 
grid  square  designator  (LTRS)  for  this  point.  The  8  digit  UTH 
grid  coordinate  for  this  point  is  then  found  and  stored  in  NBRS. 
EBCD  is  then  cailed  to  convert  the  8  digit  grid  coordinate  into 
an  alphanumeric  code  suitable  for  printing.  EBCO  expands  the  8 
digit  grid  coordinates  into  two  k  digit  variables  NBRBUF  (2) 
and  NBRBUF  (3).  These  grid  coordinates  are  then  stored  in  the 
UTH  data  table  in  locations  NBRBCD  (I)  and  NBRBCD  (2),  respectively, 
and  control  returned  to  the  calling  program.  If  an  error  is  detec¬ 
ted  during  an  attempt  to  convert  to  UTH  coordinates  an  error 
RETURN  A1  will  be  executed. 

Entry  CVNB  converts  the  8  digit  grid  coordinate  NBRS  to  an 
alphanumeric  code  suitable  for  printing  using  the  same  code  as 
was  used  by  entry  C1TA. 


P.*I.CMRT-F 


mjpjooj  p|j6  win 


BRADDOCK.  DUNN  AND  McDONALD.INC. 


1.  NAME:  CVNBM 

2.  TYPE  OF  PROGRAM:  SUBROUTINE,  LOGICAL 

3.  LANGUAGE  USED:  FORTRAN  EXTENDED 

k.  PURPOSE:  Control  the  conversion  of  a  pair  of  UTM  grid  coordinates 
to  alphanumeric  code  suitable  for  printing. 

5.  ASSUMPTIONS  AND  LIMITATIONS:  The  arguments  Nl  and-N2  are  both  8 
digit  gr • d  coordinates  in  integer  format. 

6.  ERROR  RETURNS:  None 

7.  LINKAGE  STATEMENTS  AND  DESCRIPTION  OF  ARGUMENTS: 

CALL  CVNBM  (Nl,  N2) 

Nl  -  Integer,  8  digit  UTM  grid  coordinate  of  the  lower  left  corner 
of  an  area  map 

N2  -  Integer,  8  digit  UTM  grid  coordinate  of  the  upper  right  corner 
an  area  map. 

8.  PROGRAMS  CALLING  THIS  PROGRAM:  CML,  DMHD 

9.  COMMON  VARIABLES  USED: 

A.  ZZCV3  -  NCL.  NCU 

10.  COMMON  VARIABLES  TO  BE  SET:  None 

11.  COMMON  VARIABLES  CHANGED: 

A.  ZZCV3  -  NCL,  NCU 

12;  PROGRAMS  USEO  AND  DESCRIPTION  OF  LINKAGES: 

A.  SUBROUTINES 

(1)  CALL  EBDC  (NBRS,  IDUM,  NBRBUF)  Converts  the  8  digit  integer 
UTM  grid  coordinate  into  a  two  k  digit  UTM  grid  coordinate 
representation. 
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N8RS  -  Integer,  8  digit  UTM  grid  coordinate  of  either  the 
lower  left  corner  or  upper  right  corner  of  an  area 
map . 

1 DUM  -  integer,  dummy  variable  not  used. 

NBRBUF  -  Integer,  array  of  dimension  3,  NBRBUF  (2)  being  the 
-I  digit  UTM  easting  coordinate  and  NBRBUF  (3)  being 
the  4  digit  UTM  northing  coordinate. 

13.  NOTES  ON  METHODOLOGY:  Upon  entry  to  CVNBM  a  call  is  made  to  EBCD. 

This  call  converts  the  8  digit  integer  UTM  grid  coordinate  to  a  two 
4  digit  code.  The  arguments  calculated  via  this  call  are:  IA  (2) 
integer,  4  digit  UTM  easting  grid  coordinate  for  the  lower  left  corner 
IA  (3),  integer  4  digit  UTM  northing  grid  coordinate  for  the  lower 
left  corner.  Next  NCL  (1)  is  set  to  IA  (2)  and  NCL  (2)  is  set  to  IA 
(3).  A  call  is  then  made  to  EBCD  to  convert  the  8  digit  code  N2  for 
the  upper  left  corner.  The  arguments  calculated  are  then  set  to 
NCU  (1)  and  NCU  (2),  respectively.  Control  is  then  returned  to  the 
calling  routine. 
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1 .  NAME :  EBCD 

2.  TYPE  OF  PROGRAM:  SUBROUTINE,  DATA 

3.  LANGUAGE:  FORTRAN  EXTENDED 

4.  PURPOSE:  Convert  the  8  digit  UTM  grid  coordinate  into  an  alpha¬ 
numeric  code  suitable  for  printing. 

5.  ASSUMPTIONS  AND  LIMITATIONS:  None 

6.  ERROR  RETURNS:  None 

7.  LINKAGE  STATEMENT  AND  DESCRIPTION  OF  ARGUMENTS: 

CALL  EBCD  (NBRS,  I DUM,.  NBRBUF) 

NBRS  -  8  digit  numerical  part  of  the  UTM  coordinate 
1 DUM  -  Dummy  variable. 

NBRBUF  -  integer,  array  of  dimension  3,  NBRBUF  (2)  containing  the 
easting  coordinate  and  NBRBUF  (3)  containing  the  northing 
coordinate. 

8.  PROGRAMS  CALLING  THIS  PROGRAM:  C00RT,  LABEL,  MAIN 

9.  COMMON  VARIABLES  USED:  None 

10.  COMMON  VARIABLES  TO  BE  SET:  None 

11.  COMMON  VARIABLES  CHANGED:  None 

12.  PROGRAMS  USED  AND  DESCRIPTION  OF  LINKAGES:  None 

13*  NOTES  ON  METHODOLOGY:  EBCD  encodes  and  separates  the  8  digit 
UTM  grid  coordinate  into  two  4  digit  Variables  representing  the 
easting  and  northing  of  the  8  digit  grid  coordinate.  First,  ENCODE 
is  used  to  encode  the  8  digit  integer  coordinate  and  store  this 
number  in  display  code  in  NBRBUF  (3).  NBRBUF  (3)  then  contains 
the  10  character  display  code  (2  blanks,  8  digits)  representing 
the  grid  coordinate.  This  code  is  then  shifted  right  4  charac¬ 
ters  to  drop  off  the  last  4  characters  of  the  word,  and  then 
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shifted  right  6  characters  and  stored  in  NBRBUF  (2).  This  pro¬ 
cedure  stores  a  10  digit  display  code  (k  characters,  6  blanks) 
representing  the  easting  grid  coordinate  in  NBRBUF  (2).  The  code 
in  NBRBUF  (3)  i5  then  shifted  left  6  characters  and  stored  in 
NBRBUF  (3).  This  procedure  stores  a  10  digit  display  code  (4 
character,  6  blanks)  representing  the  northing  grid  coordinate 
in  NBRBUF  (3).  Control  then  returns  to  the  calling  program. 
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1.  NAME:  ELEV 

2.  TYPE  OF  PROGRAM:  SUBROUTINE,  MATHEMATICAL 

3.  LANGUAGE:  FORTRAN  EXTENDED 

4.  PURPOSE:  Perform  s  bilinear  interpolation  to  determine  terrain 
elevation. 

5.  ASSUMPTIONS  AND  LIMITATIONS:  None 

6.  ERROR  RETURNS:  A  RETURN  A1  will  be  executed  if  |JNEfT  is  unable 
to  access  terrain  data  for  the  point  under  consideration. 

7.  LINKAGE  STATEMENT  AND  DESCRIPTION  OF  ARGUMENTS: 

CALL  ELEV  (XS,YS,Z),  RETURNS  (Al) 

XS  -  Relative  easting  of  the  point  under  consideration. 

YS  -  Relative  northing  of  the  point  under  consideration. 

Z  -  MSL  altitude  of  the  point  under  consideration. 

Al  -  Exit  taken  if  the  condition  listed  under  6  (above)  occurs. 

8.  PROGRAMS  CALLING  THIS  PROGRAM:  MAIN 

9.  COMMON  VARIABLES  USED: 

A.  ZZCV2  -  LTRS,  NBRS 

10.  COMMON  VARIABLES  TO  BE  SET: 

A.  ZZCV2  -  LTRS,  NBRS 

11.  COMMON  VARIABLES  CHANGED:  None 

12.  PROGRAMS  USED  AND  DESCRIPTION  OF  LINKAGES: 

A.  SUBROUTINES 

(1)  CALL  0NEPT(XS,YS,H)  ,  RETURNS (Al)  -  Calculate  the  terrain 
elevation  at  the  specified  grid. point. 

XS  -  Relative  easting  of  the  grid  point  under  consideration. 
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YS  -  Relative  northing  of  the  grid  point  under  consideration. 

H  -  MSL  altitude  of  the  terrain  at  the  grid  point  under 
consideration. 

A1  -  Control  is  transferred  to  this  statement  number  if 
an  error  is  detected  in  0NEPT. 

13-  NOTES  ON  METHODOLOGY:  Upon  entry  to  ELEV,  the  UTM  coordinates  of 
the  point  under  consideration  are  stored  in  local  variables  (LTRSV, 
NBRSV)  and  the  lateral  data  rate  is  set  to  500  meters.  The 
easting  and  northing  of  the  closest  western  and  southern  grid  lines 
are  calculated.  The  easting  (XREM)  and;  northing  (YREM)  from  these 
grid  lines  to  the  point  of  interest  and  the  relative  easting. (XP) 
and  northing  (YP)  of  the  intersection  of  these  grid  lines  are 
calculated.  0NEPT  is  then  called  four  times  to  calculate  the 
MSL  altitude  at  each  of  the  four  corners  of  the  square  formed  by 
the  intersection  of  the  four  grid  lines  closest  to  the  point  of 
interest.  These  altitudes  are  stored  in  local  variables  HI,  H2,  H3, 
and  H4.  The  terrain  elevation  (z)  is  then  calculated  by  bilinear 
interpolation  and  control  returned  to  the  calling  program. 
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1.  NAME:  ERR 

2.  TYPE  OF  PROGRAM:  SUBROUTINE,  OUTPUT 

3.  LANGUAGE:  FORTRAN  EXTENDED 

i».  PURPOSE:  Output  to  the  printer  an  error  message  Indicating  the 

subroutine  and  statement  number  where  an  error  occurred.  Also  Issues 
CALL  PDUHP  to  dump  variables  In  this  routine. 

5.  ASSUMPTIONS  AND  LIMITATIONS:  None 

6.  ERROR  RETURNS:  None 

7.  LINKAGE  STATEMENT  AND  DESCRIPTION  OF  ARGUMENTS: 

CALL  ERR  (NAME,  LAST,  NSTAT) 

NAME  -  The  nane  of  the  subroutine  where  error  was  detected. 

LAST  -  End  address  of  the  storage  area  to  be  dumped  from  subroutine 
NAME. 

NSTAT  -  Statement  number  In  subroutine  NAME  where  error  occurred. 

8.  PROGRAMS  CALLING  THIS  PROGRAM: 

ABSCOR  MAIN 

LABEL  PMAP 

9.  COMMON  VARIABLES  USED:  None 

10.  COMMON  VARIABLES  SET:  None 

11.  COMMON  VARIABLES  CHANGED:  None 

12.  PROGRAMS  USED  AND  DESCRIPTION  OF  LINKAGES:  None 

13.  NOTES  ON  METHODOLOGY:  Subroutine  ERR  prints  the  name  of  the 
subroutine  and  statement  number  where  an  error  occurred.  Then  a  call 
Is  made  to  PDUHP  for  a  dump  of  the  variables  used  by  this  subroutine. 
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1.  NAME:  JTRK 

2.  TYPE  OF  PROGRAM:  FUNCTION,  LOGICAL 

3.  LANGUAGE:  FORTRAN  EXTENDED 

4.  PURPOSE:  Find  the  checkpoint  code  from  the  track  data  array. 

5.  ASSUMPTIONS  AND  LIMITATIONS:  None 

6.  ERROR  RETURNS:  None 

7.  LINKAGE  STATEMENT  AND  DESCRIPTION  OF  ARGUMENTS: 

X  -  JTRK  (NTK, I ) 

NTK  -  Integer,  track  index  of  the  point  of  interest. 

I  -  Integer,  leg  index  of  the  point  of  interest. 

8.  PROGRAMS  CALLING  THIS  PROGRAM:  None 

9.  COMMON  VARIABLES  USED:  None 

10.  COMMON  VARIABLES  TO  BE  SET:  None 

11.  COMMON  VARIABLES  CHANGED:  Hone 

12.  PROGRAMS  USED  AND  DESCRIPTION  OF  LINKAGES: 

A.  SUBROUTINES 

(1)  CALL  SETBPT  (NTK,  I,  IW0RD) 

NTK  -  Index  of  desired  track 
I  -  Index  of  desired  point 

IW0RD  -  Block  pointer  address  for  the  desired  point 

13.  NOTES  ON  METHODOLOGY:  Upon  entry  to  JTRK,  the  block  pointer  word 
for  this  track  point  is  calculated  by  calling  SETBPT.  JTRK  is 
then  set  to  the  checkpoint  code  for  this  track  point.  The  debug 
flag  (DB)  is  then  tested.  If  the  debug  flag  is  TRUE,  PDUMP  is 
called  to  dump  the  track  data  array  and  control  returns  to  the 
calling  program. 
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j  ».  NAME:  LABEL 

/ 

V  . 

2.  TYPE  OF  PROGRAM:  SUBROUTINE,  LOGICAL 
!  3.  LANGUAGE:  FORTRAN  EXTENDED 

k.  PURPOSE:  Determine  the  UTM  grid  square  letter  for  PMAP  label. 

\ 

*  5.  ASSUMPTIONS  AND  LIMITATIONS:  None 

!  6.  ERROR  RETURNS:  A  cal!  to  ERR  and  a  RETURN  A1  will  be  executed  if  the 

following  error  is  detected: 

•  A.  CNL  cannot  find  a  legal  grid  letter  corresponding  to  X. 

7.  LINKACE  STATEMENT  AND  DESCRIPTION  OF  ARGUMENTS 

!  CALL  LABEL  (M,  X,  X0RG,  ID),  RETURNS  (Al). 

M  -  Character,  contains  the  one  letter  three  digit  grid  label 
for  the  line  X. 

X  -  Relative  easting  (ID=E)  or  northing  (ID*N)  for  which  the  UTM 
grid  letter  is  desired. 

XdRG  -  Displacement  (meters)  of  thr  origin  from  the  UTM  reference 
point  (AA) . 

ID  -  Character,  E  indicates  easting  of  line  is  desired.  N 
indicates  northing  of  line  is  desired. 

Al  -  Exit  taken  if  one  of  the  conditions  listed  under  6  (above) 
occurs . 

8.  PROGRAMS  CALLING  THIS  PROGRAM:  ABSC0R,  HAP 

9.  COMMON  VARIABLES  USED:  None 

10.  COMMON  VARIABLES  TO  BE  SET:  None 

11.  COMMON  VARIABLES  CHANGED:  None 
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12.  PROGRAMS  USED  AND  DESCRIPTION  OF  LINKAGES : 

A.  SUBROUTINES 

(1)  CALL  CLN  (LL,  N,  D) ,  RETURNS  (Al)  -  Calculates  th?. 

number  of  100  km  grid  squares  between  the  reference  point 
AA  and  specified  grid  square  designator. 


S? 


LL  -  Two  letter  UTM  grid  square  designator. 

N  -  Integer,  number  of  100  km  grid  squares  between  the 
specified  point  and  the  reference  point  (AA) . 

D  -  Indicates  whether  the  squares  to  be  counted  are 
easting  (E)  or  northing  (N) . 

Al  -  Control  is  transferred  to  this  statement  when  an 
illegal  grid  square  designator  is  detected  in  CLN. 

(2)  CALL  EBCD  (NBRS,  IDUM,  NBRBUF)  -  Converts  the  8  digit  UTM 
grid  coordinate  into  an  alphanumeric  code  suitable  for 
printing. 

NBRS  -  8  digit  UTM  grid  coordinate. 

IDUM  -  Dummy  variables  not  used. 

NBRBUF  -Array  of  dimension  3>  NBRBUF  (2)  containing  the 
coordinate  corresponding  to  easting  and  NBRBUF(3)  containing 
the  coordinate  corresponding  to  northing. 

(3)  CALL  ERR  (  NAME,  LAST,  NSTAT)  -  Prints  an  error  message 
indicating  name  of  subroutine  and  statement  number  where 
error  occurred. 

NAME  -  Name  of  subroutine  where  error  was  detected. 

LAST  -  End  address  of  the  storage  area  to  be  dumped  from 
subroutine  NAME. 

NSTAT  -  Statement  number  in  subroutine  NAME  where  error 
occurred. 
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(4)  CALL  MVC  (ICNT,  INPTR,  Z INADD,  0IPTR,  ZI0ADD)  -  Move  a 

character  string  from  the  word  Z INADD  into  the  word 

ZI0ADD 

ICNT  -  Integer,  number  of  characters  to  be  moved  from 
Z INADD  to  ZI0ADD. 

INPTR  -  Integer,  index  of  the  first  character  to  be  moved 
from  Z INADD. 

ZINADD  -  Character,  2  word  array  (20  character)  from  which 
characters  will  be  moved. 

I0PTR  -  Integer,  index  of  the  location  in  ZI0ADD  where 
the  first  character  will  be  moved. 

ZI0ADD  -  Character,  2  word  array  (20  characters)  where  the 
characters  will  be  moved. 

13.  NOTES  ON  METHODOLOGY:  Upon  entry  to  LABEL,  the  distance  from 
the  UTM  reference  point  (AA)  to  the  line  X  is  determined  for 
either  a  north-south  line  (ID=E)  for  an  east-west  line  (ID=N). 

The  number  (M)  of  100  km  squares  between  the  reference  and  this 
line  is  then  determined.  At  statement  1001,  CNL  is  called  to 
calculate  M,  the  grid  letter  corresponding  to  M.  If  this  operation 
does  not  result  in  an  error,  control  transfers  to  statement 
1002,  where  IX  is  determined  as  the  number  of  100  meter  steps  between 
AA  and  the  line  X  plus  1000.  EBCD  and  MVC  are  then  called  to  move 
the  three  digit  code  corresponding  to  either  the  easting  or  northing 
grid  coordinate  for  this  line.  Control  then  returns  to  the  calling 
program.  If  an  error  is  detected  by  CNL,  ERR  is  called  and  a  RETURN  A1 
executed. 
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1.  NAME:  MAIN 

2.  TYPE  OF  PROGRAM:  SUBROUTINE,  LOGICAL 

3.  LANGUAGE:  FORTRAN  EXTENDED 

4.  PURPOSE:  Control  the  logical  checking  and  plotting  of  FRAG1  type 


input  data. 

5.  ASSUMPTIONS  AND  LIMITATIONS: 

6. 


None 


ERROR  RETURNS:  A  call  to  ERR  and  a  RETURN  Al  will  be  executed  if 
any  of  the  following  errors  are  detected: 

A.  1100  -  An  illegal  grid  square  designator  has  been  specified 
for  the  coordinate  system  origin. 

4015  ~  An  1/0  parity  error  was  detected  in  an  attempt  to  read 
the  second  copy  of  the  TRACK  type  card  from  file  8. 

4016  -  An  end  of  file  was  encountered  in  an  attempt  to  read 
the  second  copy  of  the  TRACK  type  card  from  file  8. 

5001  -  An  1/0  parity  error  was  detected  in  an  attempt  to  read 
the  second  copy  of  the  P0INT  type  card  from  file  8. 

5002  -  An  end  of  file  was  encountered  in  an  attempt  to  read 
the  second  copy  of  the  TRACK  type  card  from  file  8. 

1250  -  An  1/0  parity  error  was  detected  in  an  attempt  to  read 
the  dispatch  card  from  file  8. 

3000  -  An  out  of  sequence  card  or  end  of  file  was  encountered 
in  an  attempt  to  read  the  second  copy  of  the  SITE  type 
card  from  file  8. 

3010  -  An  1/0  parity  error  was  detected  in  an  attempt  to  read 
the  second  copy  of  the  SITE  card  from  file  8. 

3080  -  The  number  of  sites  entered  exceeds  the  maximum  allow¬ 
able  (255). 

A  conversion  to  UTM  or  a  coding  of  grid  coordinates  has  resulted 
in  an  error  or  track  storage  has  been  exceeded. 


B. 


C. 


D. 


E. 


F. 


G. 


H. 


J. 
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7.  LINKAGE  STATEMENT  AND  DESCRIPTION  OF  ARGUMENTS: 

CALL  MAIN,  RETURNS  (Al) 

Al  -  Exit  taken  if  one  of  the  conditions  listed  under  6  (above) 
occurs. 

3.  PROGRAMS  CALLING  THIS  PROGRAM:  PMAP 

9.  COMMON  VARIABLES  USED: 

A.  ZZSITE  -  AZ1 A,  AZ2A,  IDFU,  IDFUTY,  XFU,  YFU,  ZFU 

B.  ZZPL0T  -  ID,  X,  Y 

C.  ZZTRK1  -  IDTHV,  IDTRK,  MNP 

D.  ZZFASC  -  KGZSV,  LTDA,  HU DA 

E.  ZZC0RD  -  LNBUF 

F.  ZZCV2  -  LTRS ,  NBRBCO,  NBRS 

G.  ZZMAPP  -  NLP  I  ,  SCALE 

10.  COMMON  VARIABLES  TO  BE  SET: 

A.  ZZFASC-KGZSV,  LTDA,  N0DA ' 

B.  ZZC0RD  -  LNBUF 

C.  ZZCV2  -  NBRBCD 

11.  COMMON  VARIABLES  CHANGED: 

A.  ZZSITE  -  AX1A,  AZ2A,  IDFU,  IDFUTY,  XFU,  YFU,  ZFU 

B.  ZZPL0T  -  ID,  X,  Y 

C.  ZZTRK1  -  IDTHV,  IDTRK,  MNP 

D.  ZZC0RD  -  LNBUF 

E.  ZZCV2  -  LTRS,  NBRS 

F.  ZZMAPP  -  NLP  I ,  SCALE 

12.  PROGRAMS  USED  AND  DESCRIPTION  OF  LINKAGES: 

A.  SUBROUTINES 

(1)  CALL  CVNB  (X,Y,L),  RETURNS  (Al)  (entry  point  to  C00RT) 
Converts  the  8  digit  UTM  grid  coordinates  into  display 
code  for  printing. 
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X  -  Dummy  variable  not  used. 

Y  -  Dummy  variable  not  used. 

L  -  Dummy  variable  not  used. 

Al  -  Control  is  transferred  to  this  statement  number  if 
an  error  is  detected  in  CVNB. 

(2)  CALL  C1TA  (X,Y,L),  RETURNS  (Al)  (entry  point  to  C00RT)- 
Converts  the  relative  easting  and  northing  of  the  point 
under  consideration  to  UTM  coordinates  (2  letters,  8  digits). 

X  -  Relative  easting  of  the  point  under  consideration. 

Y  -  Relative  northing  of  the  point  under  consideration.  i 

L  -  Dummy  variable  not  used  for  this  entry  point. 

Al  -  Control  is  transferred  to  this  statement  number  if 

an  error  is  detected  in  CITA.  i 

(3)  CALL  C1TN  (X.Y.L),  RETURNS  (Al)  (entry  point  to  C00RT)- 
Converts  UTM  coordinates  of  the  point  under  consideration 

to  relative  easting  and  northing.  I 

X  -  Relative  easting  of  the  point  under  consideration.  < 

| 

Y  -  Relative  northing  of  the  point  under  consideration.  • 

L  -  Dummy  variable  net  used.  * 

Al  -  Control  is  transferred  to  this  statement  number  i f  an  j 

error  is  detected  in  ClTN.  j 

i 

W  CALL  EBCD  (NBRS ,  I DUH,  NBRBUF)  -  Converts  the  8  digit  UTM  | 

grid  coordinate  into  an  alphanumeric  code  suitable  for  ; 

printing. 


I 

J- 

j 


V 


( 

i 


0£h 


NBRS  -  8  digit  UTM  grid  coordinate. 

IDUM  -  Dummy  variables  not  used. 

NBRBUF  -  Array  of  dimension  3,  NBRBUF  (2)  containing  the 
coordinate  corresponding  to  easting  and  NBRBUF  (2)  contain 
ing  the  coordinate  corresponding  to  northing. 
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(5)  CALL  ELEV  (XS,YS,Z),  RETURNS  (Al)  -  Determine  the 
terrain  elevation  at  a  specified  point. 

XS  -  Relative  easting  of  the  point  under  consideration. 
YS  -  Relative  northing  of  the  point  under  consideration. 

Z  -  Terrain  elevation  at  the  point  under  consideration. 
Al  -  Control  is  transferred  to  this  statement  number  if 
the  specified  point  lies  off  available  terrain. 

(6)  CALL  ERR  (NAME,  LAST,  NSTAT)  -  Prints  an  error  message 
indicating  name  of  subroutine  and  statement  number  where 
error  occurred. 

NAME  -  Name  of  subroutine  where  error  was  detected. 

LAST  -  End  address  of  the  storage  area  to  be  dumped  from 
subroutine  NAME. 

NSTAT  -  Statement  number  in  subroutine  NAME  where  error 
occurred. 

(7)  CALL  MAP  (Y,  X,  NY,  MAXNY ,  NCH) ,  RETURNS  (Al)  -  Prints  a 
map  containing  the  points  of  interest. 


Y  -  Array,  relative  northing  of  the  points  of  interest. 

X  -  Array,  relative  easting  of  the  points  of  interest. 

NY  -  Number  of  points  to  be  printed  on  the  map. 

MAXNY  -  Maximum  number  of  points  to  be  printed. 

NCH  -  Maximum  number  of  characters  in  site  identifier. 

Al  -  Control  is  transferred  to  this  statement  number  if  an 
error  is  detected  in  MAP. 


(8)  CALL  MVC  ( I CNT,  INPTR,  ZINADD,  I0PTR,  ZI0ADD)  -  Move  a 

character  string  from  the  word  ZINADD  into  the  word  ZI0ADD. 

I CNT  -  Integer,  number  of  characters  to  be  moved  from  ZINADD 
to  ZI0ADD. 

INPTR  -  Integer,  index  of  the  first  character  to  be  moved 
from  ZINADD. 
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ZINADD  -  Character,  2  word  array  (20  character)  from 
which  characters  will  be  moved. 

I0PTR  -  Integer,  index  of  the  location  in  ZI0ADD  where 
the  first  character  will  be  moved. 

ZI0ADD  -  Character,  2  word  array  (20  characters)  where 
the  character  will  be  moved. 

(9)  CALL  PACK  (I,  IS1TE,  ZWENG)  -  Prints  the  track  site 
engagement  data  on  file  9. 

I  -  Integer,  track  index. 

ISITE  -  Integer,  site  index. 

ZWENG  -  Logical,  site- track  engageabi I i ty  flag. 

(10)  CALL  SEARCH  (ISITE)  -  Prints  track  identifiers  for 
engageable  site-track  combinations. 

ISITE-  Integer,  site  index  for  site  Under  consideration. 

(11)  CALL  SET0RG  (X,Y,L),  RETURNS  (A1 )  (entry  point  to  C00RT) 

Set  coordinate  system  origin  to  grid  square  specified. 

X  -  Dummy  variable  not  used. 

Y  -  Dummy  variable  not  used. 

L  -  Two  letter  grid  square  designator  for  the  origin  of 
the  area  where  the  game  is  to  be  played. 

Al  -  Control  is  transferred  to  this  statement  number  if  an 
error  is  detected  in  SET0RG. 

(12)  CALL  S0RTR  (X,Y,NX,MASNX,NY,M0DE)-  Perform  sorts  on  the 
speci fied  arrays. 

X  -  Array  for  which  an  ascending  sort,  descending  sort, or 
maximum  and  minimum  of  the  array  entries  will  be 
calculated. 

Y  -  Array  whose  entries  will  be  moved  in  an  order  correspond 

ing  to  the  sort  performed  on  X. 
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NX  -  Integer,  number  of  entries  from  the  array  X  to  be 
avai lable  for  sort. 

MAXNX  *•  Integer,  maximum  number  of  entries  upon  which  a  sort 
is  to  be  performed. 

NY  -  Integer,  number  of  entries  in  array  Y  which  will  be 
moved  for  each  corresponding  move  in  the  X  array. 

M0DE  -  Integer,  mode  of  operation  for  this  routine,  1  implies 
descending  sort, 2  implies  ascending  sort,  and  3  implies 
maximum  and  minimum  entries  only  are  desired. 

(13)  CALL  STPNT  (X,  Y,  Z,  T,  J),  RETURNS  (El)  -  Stores  track 

point  data  array. 

X  -  Relative  easting  of  the  track  point  under  considera¬ 
tion. 

Y  -  Relative  northing  of  the  track  point  under  consideration. 

Z  -  Altitude  of  the  track  point  under  consideration. 

T  -  Time  required  to  reach  the  point  under  consideration. 

J  -  Checkpoint  code  for  the  point  under  consideration. 

El  -  Control  is  returned  to  this  statement  number  if  the 
track  point  data  array  overflows. 

(14)  CALL  STRK  (DUM,  DUM,  DUM,  DUM,  DUM) ,  RETURNS  (El) 

(entry  point  to  STPNT)  -  Initializes  the  block  pointer 
for  storage  of  this  tracks  points. 

DUM  -  Dummy  variable  not  used. 

El  -  control  is  returned  to  this  statement  number  if  the 
track  point  data  array  overflows. 

(15)  CALL  TT0RNR  (RAGIN,  TIZZ,  TIME,  HF,  NTK)  -  Calculates 
the  track  entry  time  into  the  specified  site  fire 
volume. 
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RAG  IN  -  Range  of  Interest  for  the  site  under  consideration. 
TIZZ  -  Earliest  time  to  be  considered  for  entry. 

TIME  -  Time  a  specified  track  enters  a  specified 
fire  volume. 

NF  -  Integer,  index  of  the  site  under  consideration. 

NTK  -  Integer,  index  of  the  track  under  cors ideration. 

B.  FUNCTIONS 

(1)  X.”  TTRK(NTK,I)  -  Finds  the  time  associated  wi th  a 
specified  track  point. 

NTK  -  Integer,  index  of  the  desired  track. 

I  -  Integer,  index  of  the  desired  point  on  the  track. 

(2)  X  =  XTRK(NTK, I )  (entry  point  to  TTRK)  -  Finds  the 
relative  eastinq  of  the  specified  track  point. 

NTK  -  Integer,  index  of  the  desired  track. 

1  -  Integer,  index  of  the  desired  point  on  the  track. 

(3)  21“  YTRK(NTK, I)  (entry  point  to  TTRK)  -  Finds  the  relative 
northing  of  the  specified  track  point. 

NTK  -  Integer,  index  of  the  desired  track. 

I  -  Integer,  index  of  the  desired  point  on  the  track. 

(4)  X  =  ZTRK(NTK, I ) (entry  point  to  TTRK)-  Finds  the  terrain 
elevation  of  the  specified  track  point. 

NTK  -  Integer,  index  of  the  desired  track. 

I  -  Integer,  index  of  the  desired  point  on  the.  track. 
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13.  NOTES  ON  METHODOLOGY:  Upon  entry  to  subroutine  MAIN  the  default 
values  for  the  plot  scale  (SCALE),  the  loyical  option  flags,  and 
UTM  coordinates  for  the  oriyin  are  initialized.  The  mathematical 
origin  is  then  calculated  by  a  call  to  subroutine  SET0RG.  Following 
initialization  of  the  default  values,  the  first  input  card  is  read 
from  file  5  and  stored  in  the  record  buffer  (IBUF).  If  this 
record  is  not  an  end-of-file  marker,  control  transfers  to  state¬ 
ment  1150  and  two  copies  of  the  above  mentioned  record  are 
punched  on  file  8.  Control  then  transfers  to  statement  1100  for 
processing  of  the  next  input  card.  When  an  end  of  file  is 
encountered  on  file  5,  control  transfers  to  statement  1200  where 
an  end  of  file  and  rewind  are  executed  for  file  8.  Next  the  site, 
system,  track,  and  error  counters  are  initialized  to  zero  and 
Che  logical  processing  flags  are  initialized. 

Statements  1250  through  1400  read  the  first  copy  of  each 
card  type  from  the  secondary  input  device,  determine  the  card's 
dispatch  number,  and  control  branching  for  the  appropriate  card 
types.  At  statement  1250,  the  card  counter  (N)  is  incremented 
and  the  first  copy  of  the  input  card  is  read  from  file  8  and 
stored  in  the  buffer  (IBUF).  If  an  end  of  file  is  encountered, 
control  transfers  to  statement  1975  for  final  data  processing. 
Otherwise,  an  \/t  parity  error  test  is  performed  and  if  an  error 
is  detected,  control  transfers  to  statement  9010  for  error  pro¬ 
cessing.  A  DO  loop  is  then  used  to  find  the  cards  dispatchinq  index. 
If  this  is  a  nondispatchable  card  type,  control  transfers  to 
statement  1300  where  an  error  message  is  printed  indicating  an 
illegal  card  type  has  been  found,  the  second  copy  of  the 
nondispatchable  card  type  is  read  from  file  8,  the  error  counter 
is  incremented,  and  control  transfers  to  statement  1250  where 
the  next  card  type  is  processed.  Dispatchable  indices  of  I,  2, 

3,  4,  and  5  corresponding  to  card  types  of  SWIT,  SITE,  TRAC, 

P0IN,  and  SYST,  cause  a  corresponding  transfer  of  control  to 
statements  2000,  3000,  4000,  5000,  and  6000. 

P.4.HAIN-H  _ 


71< 


BRADDOCK.  DUNN  AND  McDONALD. INC. 


At  statement  1975,  the  logical  flag  ZSILL  is  tested  to 
determine  If  any  SITE  cards  have  been  processed.  If  ZSILL  is 
TRUE,  no  site  cards  have  been  processed  and  control  transfers 
to  statement  1985-  If  ZSILL  is  FALSE  and  plotting  is  desired, 
a  plot  of  these  sites  will  be  printed  and  control  returns  to  the 
calling  program.  At  statement  1985  logical  flag  ZNOPLT  is 
tested  to  determine  if  plotting  is  desired.  If  ZN0PLT  is  FALSE, 
plotting  is  desired  and  the  appropriate  track  points  will  be 
plotted.  If  no  plotting  is  desired,  control  transfers  to 
statement  1995  where  the  no  engageabi 1 i ty  flag  is  tested.  If 
no  engageabi 1 i ty  is  to  be  played,  control  returns  to  the  calling 
program.  Otherwise,  if  system  sites  and  tracks  were  entered, 
control  transfers  to  statement  8000  where  the  engageabi 1 i ty  test 
is  performed.  Control  then  transfers  to  the  calling  program. 

At  statement  2000  the  second  copy  of  the  SWITCH  card  is 
read  from  the  secondary  input  device  (file  8).  Variables  read 
from  the  SWITCH  card  are:  KGZ,  UTM  grid  zone  designator;  K0R, 
coordinate  system  origin;  NLPI ,  number  of  lines  per  inch;  SCALE, 
plotting  scale  factor;  RNMAX,  maximum  leg  length;  ANG0,  maximum 
turn  angle;  ZTER,  logical  terrain  flag;  ZPUNCH,  logical  site 
punch  flag;  ZN0PLT,  logical  no  plot  flag;  ZN0ENG,  logical  no 
engageabi 1 i ty  flag,  and  ZN0SEC.  logical  no  sector  flag.  Following 
statement  2000,  the  cosine  of  the  maximum  turning  angle  is 
calculated,  the  number  of  lines  per  page  (NLPP)  is  calculated, 
and  the  mathematical  origin  is  determined  by  a  call  to  SET0RG. 

At  statement  2100  the  logical  terrain  flag  is  tested.  If 
terrain  is  not  to  be  considered  by  this  program  (ZTER  =  FALSE), 
control  transfers  to  statement  1250.  If  terrain  is  tc  be  con¬ 
sidered,  the  terrain  directory  is  read  from  the  direct  access 
device  (file  k)  and  control  transfers  to  statement  1250. 
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At  statement  1975,  the  logical  flag  ZSILL  is  tested  to 
determine  if  any  SITE  cards  have  been  processed.  If  ZSILL  is 
TRUE,  no  site  cards  have  been  processed  and  control  transfers 
to  statement  1935.  If  ZSILL  is  FALSE  and  plotting  is  desired, 
a  plot  of  these  sites  will  be  printed  and  control  returns  to  the 
calling  program.  At  statement  1985  logical  flag  ZNOPLT  is 
tested  to  determine  if  plotting  is  desired.  If  ZN0PLT  is  FALSE, 
plotting  is  desired  and  the  appropriate  track  points  will  be 
plotted,  if  no  plotting  is  desired,  control  transfers  to 
statement  1995  where  the  no  engageabi 1 i ty  flag  is  tested.  If 
no  engageabi 1 i ty  is  to  be  played,  control  returns  to  the  calling 
program.  Otherwise,  if  system  sites  and  tracks  were  entered, 
control  transfers  to  statement  8000  where  the  engageabi lity  test 
is  performed.  Control  then  transfers  to  the  calling  program. 

At  statement  2000  the  second  copy  of  the  SWITCH  card  is 
read  from  the  secondary  input  device  (file  8).  Variables  read 
from  the  SWITCH  card  are:  KGZ,  UTM  grid  zone  designator;  K0R, 
coordinate  system  origin;  NLPI,  number  of  lines  per  inch;  SCALE, 
plotting  scale  factor;  RNMAX,  maximum  leg  length;  ANG0,  maximum 
turn  angle;  ZTER,  logical  terrain  flag;  ZPUNCH,  logical  site 
ounch  flag;  ZNflPLT,  logical  no  plot  flag;  ZN0ENG,  logical  no 
engageabi 1 i ty  flag,  and  ZN0SEC,  logical  no  sector  flag.  Following 
statement  2000,  the  cosine  of  the  maximum  turning  angle  is 
calculated,  the  number  of  lines  per  page  (NLPP)  is  calculated, 
and  the  mathematical  origin  is  determined  by  a  call  to  SET0RG. 

At  statement  2100  the  logical  terrain  flag  is  tested,  if 
terrain  is  not  to  be  considered  by  this  program  (ZTER  =  FALSE), 
control  transfers  to  statement  1250.  If  terrain  is  to  be  con¬ 
sidered,  the  terrain  directory  is  read  from  the  direct  access 
device  (file  4)  and  control  transfers  to  statement  1250. 
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At  statement  6000,  the  processing  of  SYSTEM  type  data 
cards  begins,  first  the  system  counter  (H0SYS)  is  incremented 
and  the  second  copy  of  the  SYSTEM  type  card  is  read  from  the 
secondary  input  device  (file  8).  The  system  identifier  (lOSYS) 
and  the  maximum  range  of  interest  (RANMAX)  are  read  from  the 
SYSTEM  type  card  and  stored  in  the  system  data  table.  Control 
then  transfers  to  statement  1250. 

At  statement  3000,  the  processing  of  SITE  type  data  cards 
begins.  If  this  is  the  first  SITE  type  card  to  be  processed, 
the  site  counter  (NFU)  is  initialized  to  zero.  A  test  is  then 
performed  on  logical  flag  ZSILL  to  determine  if  this  SITE  card 
is  being  read  out  of  sequence.  If  this  SITE  card  is  being  read 
out  of  sequence,  control  transfers  to  statement  9200  for  error 
processing.  If  this  card  is  not  out  of  sequence  the  logical 
flags  ZERCVT,  ZSP,  and  Z0N  are  then  initialized.  Next,  the  line 
buffer  is  initialized  and  the  site  altitude  is  set  to  zero.  At 
statement  30i0,  the  second  copy  of  the  SITE  type  data  card  is 
read  from  file  8.  Variables  contained  on  the  SITE  type  data 
card  are:  IFTYP,  type  of  system;  IDF,  site  identifier;  ITRS, 
two  letter  grid  square  designator;  HBRS,  eight  digit  UTM  grid 
coordinate;  ZF,  altitude  of  the  site  above  MSI;  AS1,  the  left 
,  edge  of  sector;  and  AZ2,  the  right  edge  of  sector.  If  an  end 
of  file  is  encountered  or  an  1/0  parity  error  is  detected, 
control  transfers  to  statements  9200  and  9220,  respectively, 
for  error  processing.  If  no  1/0  errors  are  detected, control 
transfers  to  statement  3015  where  the  relative  easting  and 
northing  of  the  site  locations  are  determined  by  call  to 
C1TN.  If  terrain  is  being  considered,  the  terrain  elevation 
(Z)  is  calculated  by  a  call  to  ELEV  and  the  difference  between 
the  calculated  and  specified  terrain  altitudes  is  calculated. 

If  this  site  location  is  found  to  be  off  the  available  terrain 
area,  ELEV  sets  the  logical  flag  Z0W  to  FALSE.  The  sites 
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description  is  then  printed  in  the  site  description  table  along 
with  a  note  describing  whether  the  site  is  off  the  available 
terrain.  A  test  is  then  performed  and  if  the  site  altitude 
specified  is  below  the  terrain  elevation,  an  appropriate  error 
message  is  printed.  At  statement  3019  the  site  counter  (NFU)  is 
tested.  If  the  number  of  sites  is  zero  no  redundancy  test  is 
required  and  control  transfers  to  statement  3080.  Otherwise,a 
DO  loop  is  entered  which  tests  the  site  name  against  the  site 
names  already  entered  in  the  site  description  table.  If  this 
site  name  is  unique,  cor.  >■•)!  transfers  to  statement  3050.  If 
a  redundant  site  name  is  found,  control  transfers  to  statement 
3030  where  the  specified  and  redundant  site  locations  are  tested. 

If  the  site  locations  are  the  same,  a  message  is  printed  indicating 
that  a  duplicate  SITE  card  has  been  found  and  control  transfers 
to  statement  1250  where  the  next  card  type  will  be  processed. 
Following  statement  3040,  a  logical  change  flag  (ZCH)  is  set 
TRUE  and  an  error  message  printed  to  indicate  that  a  duplicate 
card  ID  has  beer,  found.  At  statement  3050,  a  DO  loop  is  used  to 
determine  if  this  site  location  is  unique.  If  this  site  location 
is  unique,  control  transfers  to  statement  3080.  Otherwise, 
control  transfers  to  statement  3070  where  an  error  message  is 
printed  indicating  that  two  sites  have  been  found  with  the  same 
location,  the  change  flag  ZCH  is  set  to  TRUE  and  control 
transfers  to  statement  1250.  Following  statement  3080,  the  site 
description  is  entered  in  the  site  description  table  and  control 
transfers  to  statement  1250. 

TRACK  type  data  cards  are  processed  starting  at  statement 
4000.  At  statement  number  4000,  the  logical  flag  ZSITE  is 
tested  to  determine  if  site  cards  have  been  previously  processed. 

If  site  cards  have  been  previously  processed,  the  plotting  flag 
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(ZN0P10T)  is  tested  to  determine  if  plotting  is  desired.  If  no 
sites  have  been  processed  or  no  plotting  is  desired,-  control 
transfers  to  statement  4005.  Otherwise,  the  sites  are  sorted  and 
subroutine  HAP  is  cal  led  to  plot  the  desired  site  map.  Fol  lowing 
statement  4005,  logical  flags  ZS1T  and  ZSILL  are  reset  to 
disallow  further  site  processing  and  track  counter  (NTK)  is 
incremented.  At  statement  4015,  the  second  copy  of  the  track 
type  data  card  is  read  from  file  8.  if  an  end  of  file  or  i/0 
parity  error  is  detected,  control  transfers  to  statement  9004 
or  9006,  respectively,  for  the  appropriate  error  processing.  The 
variables  read  from  the  TRACK  type  data  card  and  stored  in  the  track 
description  tables  are:  IDTHV,  penetrator  identifier;  IDTRK, 
track  identifier;  and  MNP,  number  of  track  points.  A  DO  loop 
is  then  used  to  test  for  redundant  track  identifiers  in  the  track 
description  table.  If  a  redundant  track  identifier  is  found, 
control  transfers  to  statement  5040  where  an  error  message  is 
printed  to  indicate  that  a  redundant  track  identifier  has  been 
found.  Control  then  transfers  to  statement  4030.  At  statement 
4030  a  test  is  performed  to  determine  if  the  number  of  points 
(MN)  processed  for  the  previous  track  is  less  than  zero.  If  no 
points  were  processed,  control  transfers  to  statement  4100 
where  an  error  message  is  printed.  If  noints  were  processed  for 
the  previous  track,  STRK  is  called  to  determine  the  initial 
storage  address  for  the  data  on  this  track.  A  test  is  then 
performed  on  the  track  counter  (NTK),  and  control  is  transferred 
to  statement  4200  if  this  is  the  first  track  card  processed, 
if  this  is  not  the  first  track  to  be  processed,  the  actual 
point  count  for  this  track  and  the  logical  track  point  error 
flags  are  tested.  If  an  error  is  detected  in  either  one  of 
these  values,  an  error  message  will  be  printed.  Otherwise,  the 
no  plot  flag  is  tested  and  if  plotting  is  desired,  the  track 
points  will  be  sorted  and  a  track  map  will  be  plotted.  At 
stotement  4200,  track  data  for  this  track  will  be  printed. 
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Following  statement  4300  the  logical  error  flags  and  the  point 
cot  ter  for  this  track  will  be  initialized  and  control  then 
returns  to  statement  1250  for  processing  of  the  next  card  type. 

Processing  of  POINT  type  cards  begins  at  statement  5000. 

Following  statement  5000,  the  point  counter  (MN)  is  incremented 
and  the  track  time  (TTK)  is  initialized  to  zero.  At  statement 
5001.  the  second  copy  of  the  POINT  type  data  card  is  read  from 
file  8.  I f  an  end  of  file  or  an  1/0  parity  error  is  detected, 
control  will  transfer  to  statement  9008  or  9006,  respectively, 
for  error  processing.  The  variables  read  from  the  POINT  type 
card  are:  LTRS,  two  letter  grid  square  designator;  NBRS, 
eight  digit  UTM  grid  coordinate;  ZTK,  altitude  of  the  track  above 
main  sea  level;  VEL,  velocity  of  the  vehicle  approaching  that  point; 
IDCP,  check  point  code;  IDTK,  track  identifier;  IP,  sequence  number  of 
the  point  on  the  track;  and  I  DTP,  track  point  identifier.  After 
the  POINT  type  card  has  been  successfully  read  from  file  8,  C1TN 
is  called  to  convert  the  UTM  coordinates  of  this  point  to  relative 
coordinates  and  CVNB  is  called  to  code  the  UTM  grid  coordinates 
in  a  form  suitable  for  printing.  The  targeted  track  logical 
flags  are  then  initialized  and  a  test  is  performed  to  see  if  the 
track  point  identifier  has  been  left  blank.  If  the  track  point 
identifier  has  been  left  blank,  this  track  point  is  not  targeted 
for  any  site  and  control  transfers  to  statement  5026.  If  a  track 
point  identifier  has  been  specified,  a  loop  is  then  entered 
which  tests  each  site  identifier  against  the  track  point  identifier. 
When  it  is  determined  that  a  track  is  targeted  for  particular 
site,  control  transfers  to  statement  9310.  If  no  targeted  site 
is  found,  logical  flag  ZTTENG  is  set  to  FALSE  and  control 
transfers  to  statement  5026.  At  statement  9310  the  square  of 
the  horizontal  range  between  the  site  location  and  track  point  is 
calculated  and  tested  to  determine  if  this  track  point  is  within 
25  meters  of  the  site  location.  If  this  track  point  lies  within 
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25  meters  of  the  site  location,  logical  flag  ZWENG  is  set  to 
TRUE.  At  statement  5026,  PACK  is  called  to  punch  a  track/sit* 
targeting  record  on  file  9.  If  this  is  the  first  point  to  jk 
processed  for  this  track,  a  track  table  heading  is  printed  and 
this  leg's  slant  range  (RTK)  and  ground  track  (GRNG)  are 
initialized.  Then  control  transfers  to  statement  5010.  If  this 
is  not  the  first  point  to  be  processed  on  this  track,  GRNG,  RTK, 
and  the  time  required  for  the  vehicle  to  reach  this  point  on  the 
track  (TTK)  are  calculated  following  statement  5009.  If  required, 
a  skip  to  the  top  of  the  next  page  is  issued  and  a  new  table 
heading  is  printed  following  statement  5010.  At  statement  5011 
the  track  point  data  are  printed  in  the  track  data  table.  If 
the  site  was  not  found  with  the  specified  track  point  identifier, 
an  appropriate  message  is  then  printed  for  this  track.  If 
terrain  is  to  be  tested,  a  test  is  then  performed  to  determine 
if  this  point  lies  on  available  terrain.  If  this  point  lies 
off  the  available  terrain,  an  appropriate  error  message  is 
printed.  Statements  5013  through  5019  test  this  leg  for  four 
different  types  of  errors.  If  this  point  track  identifier 
is  found  to  differ  from  the  actual  identifier  for  this  track, 
logical  flag  ZTEKID  is  set  to  TRUE  and  an  appropriate  error 
message  is  printed.  If  the  point  sequence  number  does  not 
correspond  to  the  actual  sequence  number  for  this  point,  logical 
flag  ZETKPT  is  set  to  TRUE  and  an  appropriate  error  message  is 
printed.  If  this  leg's  length  is  less  than  10  kilometers, 
logical  flag  ZSP  is  set  to  TRUE  and  a  message  is  printed  indicating 
that  this  leg  is  unusually  short.  If  RTK  is  greater  than  the 
maximum  range  specified  on  the  SWITCH  card,  a  message  is 
printed  indicating  that  an  unusually  long  leg  has  been  found 
and  logical  flag  ZSP  is  set  to  TRUE.  If  the  number  of  points 
processed  is  greater  than  two,  a  test  is  then  performed  to 
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determine  if  the  turning  angle  defined  by  the  previous  track 
points  is  greater  than  the  maximum  expected  turning  angle.  If 
the  maximum  expected  turning  angle  is  exceeded,  ZSP  is  set  to 
TRUE  and  an  appropriate  message  is  printed.  Following  statement 
5020  this  point  is  stored  in  the  track  data  array  and  the  point 
sequence  number  is  stored  in  variable  ID.  Control  is  then 
returned  to  statement  1250.  Statements  5030  through  5900  set 
error  flags  and  print  appropriate  error  messages  if  storage, 
redundant  tracks,  or  conversion  errors  are  detected  for  this 
point.  Then  control  transfers  to  statement  1250  for  processing 
of  the  next  card  type. 

Following  statement  8000,  the  heading  for  the  site 
engageabi 1 i ty  table  is  printed.  A  major  loop  is  then  entered 
which  determines  which  tracks  are  engageable  by  which  sites. 

First  the  site  system  index  ( I J )  is  determined.  If  the  system 
type  cannot  be  found  for  this  site,  an  error  message  is  printed 
and  control  transfers  to  statement  8100.  Once  the  site  system 
type  has  been  determined,  the  site's  range  of  interest  is  stored 
in  variable  RANG.  A  DO  loop  is  then  used  to  determine  the 
tracks  which  intersect  this  fire  volume.  For  each  track,  TT0RNG 
is  called  to  determine  the  time  that  the  vehicle  fiying  this 
track  will  enter  the  fire  volume.  If  this  time  is  negative, 
the  track  does  not  enter  the  fire  volume  and  control  transfers 
to  statement  8020.  Otherwise,  the  intersection  counter  and 
intersection  flags  are  set  and  the  track  ID  is  entered  in  the 
site  engageabi 1 i ty  table.  An  appropriate  engageabi 1 i ty  message 
is  then  printed  for  the  site.  Following  statement  8100,  a  DO 
loop  is  entered  which  prints  a  table  containing  the  identifiers 
for  all  the  nonengageable  tracks.  If  all  the  tracks  are 
nonengageable  by  all  the  sites,  an  appropriate  message  is  printed 
following  statement  8200.  A  major  DO  loop  is  then  entered 
which  determines  which  tracks  are  targeted  for  which  sites  and 


P.4. MAIN-0 


73< 


BRADDOCK.  DUNN  AND  MCDONALD, INC. 


prints  a  table  of  targeted  sites.  First  the  logical  flags  ZTENG 
is  set  TRUE,  if  the  track  has  a  point  which  comes  within  25 
meters  of  the  site  location.  Next  a  test  is  performed  and  if 
the  track  point  lies  within  the  site's  range  of  interest,  the 
track  identifier  is  stored  in  the  table  ITP.KSW.  A  message  is 
then  printed  indicating  which  tracks  are  targeted  for  each  site 
and  control  is  returned  to  the  calling  program. 
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1.  NAME:  MAP 

2.  TYPE  OF  PROGRAM:  SUBROUTINE,  OUTPUT 
3;  LANGUAGE:  FORTRAN  EXTENDED 

4.  PURPOSE:  Print  'a  map  containing  the  points  of  interest. 

5.  ASSUMPTIONS  AND  LIMITATIONS:  The  coordinates  of  the  points  of 
interest  have  be.en  appropriately  loaded  into  the  arrays  in  the 
proper  form. 

6.  ERROR  RETURNS:  A  call  to  ERR  and  a  RETURN  A1  will,  be  executed 
if  either  of  the  'following  errors  are  detected : 

A.  2000-  ABSC0R  is  unable  to  generate  the  required  abscissa 

coordinate  labels. 

B.  2900  -  LABEL  is  unable  to  generate  the  ordinate  labels. 

7.  LINKAGE  STATEMENT  AND  DESCRIPTION  OF  ARGUMENTS: 

CALL  Map  (Y,  X,  NY,  MAXNY ,  NCH) ,  RETURNS  (A1): 

Y  -  Array,  relative  northing  of  the  points  of  interest. 

X  -  Array,  relative  easting  of  the  points  of  interest. 

NY  -  Number  of  points  to  be  printed  on  the  map. 
fMAXNY  -  Maximum  number  of  points  to  be  printed. 

NCH  -  Maximum  number  of  characters  in  site  identifier. 

A1  -  Exit  taken  if  one /of  the  conditions  listed  under  6(above)  occurs. 
8..  PROGRAMS  CALLING  THIS  PROGRAM:  MAIN 
9.  COMMON  VARIABLES  USED: 

A.  ZZC0RD  -  DELTA,  ESTEDG,  I0RBUF,  LNBUF  (equivalenced  to 

ZLNBUF) ,  WSTEDG 

B.  ZZMAPP  -  NLP  I ,  SCALE  -  XMAX,  XMIN,  YMAX  ,  YMIN 

C.  ZZCVI  -  Y0RG 
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10.  COMMON  VARIABLES  TO  BE  SET: 

A.  ZZC0RD  -  I0RBUF,  LNBUF  (equ 1 va 1 enced  to  ZLNBUF) 

B.  ZZMAPP  -  NLPI,  SCALE,  XMAX,  XMIN,  YMAX,  YMIN 

C.  ZZCV1  -  Y0RG 

11.  COMMON  VARIABLES  CHANGED: 

A.  ZZC0RD  -  DELTA,.  ESTEDG,  LNBUF  (equ I va fenced  to  ZLNBUF) ,  WSTEDG 

B.  ZZMAPP  -  XMAX,  YMIN 

12.  PROGRAMS.  USED  AND  DESCRIPTION  OF  LINKAGES: 

A.  SUBROUTINES 

(1) '  CALL  ABSC0R  (l START,  ISPC,  N),  RETURNS  (Al) 

(START  -  Integer,  tha  starting  position  for  the 
abscissa  coordinates  in  the  print  array. 

ISPC  -  The  number  of  blanks  between  coordinate  labels. 

N  -  Number  of  labels  to  be  printed  across  the  page. 

Al  -  Control  transfers  to  this  statement  number  if  an 
error  is  detected  in  ABSC0R. 

(2)  CALL,  ERR  (NAME,  LAST,  NSTAT)  -  Prints  an  error  message 
Indicating  name  of  subroutine  and  statement  number 
where  error  occurred. 

NAME  -  Name  of  subroutine  where  error  was  detected. 

LAST  -  End  address  of  the  storage  area  to  be  dumped  from 
subroutine  NAME. 

NSTAT  -  Statement  number  in  subroutine  NAME  where  error 
occurred i 

(3)  CALL  LABEL  (M,  X,  XORG,  ID),  RETURNS  (Al) 

M  -  Character,  contains  the  one  letter  three  digit 
grid*  label  for  the  line  X.  . 
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X  -  Relative  easting  (lD*E)  or  northing  (ID«N)  for 
which  the  UTM  grid  letter  is  desired. 

X0RG  -  Displacement  (meters)  of  the  origin  from  the  UTM 
reference  point  (AA) . 

ID  -  Character,  E  indicates  easting  of  strip  is 
desired.  N  indicates  northing  of  line  is 
desi red. 

Al  -  Control  is  transferred:  to.  this,  statement  when 
an  illegal  grid  letter  is  detected  by  CLN. 

(*).  CALL  MVC  ( I CNT ,  INPTR,  Z INADD,  I0PTR,  ZI0ADD)  -Move 
a  character  string  from  the  word  Z INADD  into  the 
word  ZI0ADD. 

I CNT  -  Integer,  number  of  characters  to  be  moved 
from  Z INADD  to  ZI0ADD. 

INPTR  -  Integer,  index  of  the  first  character  to 
be  moved  from  Z INADD. 

ZINADD  -  Character,  2  word  array  (20  character) 
from  which  characters  will  be^moved. 

I0PTR  -  Integer,  index  of  the  location  in  ZI0ADD 
where  the  first  character  will  be  moved. 

ZI0ADD  -  Character,  2  word  array  (20  characters) 
where  the  character  wi 1 1  be  moved. 

13.  NOTES  OF  METHODOLOGY:  Upon  entry,  MAP  sets  the  scale  factor 
(I SCALE),  calculates  the  easting  scale  factor  (DELTAX) , 
calculates  the  northing  scale  factor  (DELTAY),  calculates 
the  ordinate  step  size,  and  calculates  the  map  boundries. 

Next,  the  strip  counter  ( I  STRIP) ,  which  counts  the;  number 
of  northrsouth  strips  required  to  cover  the  desired  area 
map,  is  initialized  to  zero. 
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Next,  a  loop  is  entered  which  calculates  and  plots 
the  north-south  strips  required  to  coyer  the  above  mentioned 
map.  Following  statement  1000,  the  eastern  edge  (ESTEDG) , 
and  western  edge  (WSTEDG)  of  the,  strip  are  calculated; 

The,  western  edge  of  the  next  strip  (XNXT)  is  calculated  and 
the  logical  flag  (ZRLAB)  is  set  TRUE  i f  no  more  strips 
are  required  to  cover  the  desired  map.  ABSC0R  is  then  called 
to  store  the  abscissa  labels  for  this  strip  in  buffer 
I0RBUF.  the  strip  counter  (ISTRIP)  is  then  incremented 
and  tested..  If  the  strip  under  consideration  |s  the  first, 
strip  to  be  processed  for  this  map,  a  page  eject  is  executed 
and  the  map  header  data  printed.  Otherwise;  a  skip  to  the 
top  of  the  next  page  is  executed.  At  statement  2020, 
the  one  letter  three  digit  abscissa  label  is  printed  across 
the  page.  Next  the  site  counter  (!)„  the  northing  of  the 
line  to  be  printed  (YAXIS)  and  the  line  counter  (KG)  are 
initialized  for  this  strip.  At  statement  2100  the 
northern  and  southern  boundaries  of  the  area  which  this  line 
of  print  represents  are  calculated.  Next  the  JM)  loop  is  used 
to  clear  the  line  buffer.  Logical  flag  ZLABEL  then  is  set 
to  TRUE  if  an  ordinate  label  is  required  for  this  line. 

The  grid  for  this  line  is  next  created.  First  the  character 
is  placed  in  the  line  buffer  in  locations  corresponding 
to  the  abscissa  labels.  If  an  ordinate  label  is  required 
(ZLABEL  *  TRUE),  the  line  buffer  is  then  filled  to  print 
a  complete  line  of  grid  points.  At  statement  2500,  the  si'tc 
counter  is  then  tested.  If  all  the  sites  have  been  previously 
mapped,  controltransfers  to  statement  2800.  Following 
statement  2500,  logical  flag  Z0UTRG  is  set  TRUE,  if  site  I  does 
not  lie  on  the  area  represented  by  this  line  and  logical 
flag  ZC0INC  is  set  TRUE  if  the  next  site  in  the  site  table 
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(1+1)  lies  on  this  area.  If  this  point  (|)  does  not  lie 
on  this  area,  control  transfers  to  statement  2800.  Otherwise, 
a  test  is  performed*  to  determine  if  the  point  I  dies,  within 
the  strip.  If  this  point  lies  within  this  strip 'and  wi thi n 
this  line  area,  the  point  location  (L0C)  in  the  line  buffer 
is  then  calculated.  If  this  location  lies  outside  the  printer 
page  an  error  message  is  printed  and;  L0C  is  reset  to  l. 
Following  statement  2600  the  line  counter  is  incremented 
by  seven  to  provide  the  line  offset,  the  number  of  characters 
identifying  the  point  is  calculated- and  stored  in  the  line 
buffer  in  the  appropriate  location  for  this  point.  The  point 
counter'  is  then  incremented  and  logical  flag  ZC0INC  is  tested 
to  determine  if  another ;point  lies  on  this  line  area.  If 
it  is  determined’ that  another  point  is  to  be  printed  on  this 
line,  control  transfers  to  statement  2500.  A  test  is  then 
performed  to  determine  if  an  ordinate  label  is  required  for 
this  line..  If  a  label  is  fequi  red  for  this  line,  LABEL  is 
then  cal.ied- to  find; this;  axis  label  and  HVC  is  used  to  store 
this  'label  in  the  line  buffer.  This  line  is  then  printed  on 
the  map 'and  the  northing  is  indexed  to  the  next  line,  and 
control  is  transferred  to  statement  2100  if  printing,  for  this 
strip  is  not  complete.  Otherwise,  the  abscissa  labels  are 
printed  on  the  end  of  the  map  and  control  returns  to  statement 
1000  if  it  is  required  to  process  another  strip. 

Otherwise,  the  normal  return  is  executed  to  the  calling 
program. 
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1.  NAME:  MAXN  (other  entry  points:  MINN) 

2.  TYPE  OF  PROGRAM:  SUBROUTINE,  LOGICAL 

3.  LANGUAGE:  FORTRAN  EXTENDED 

4.  PURPOSE:  Determine  the  array  entry  with  the  greatest  (smallest)  value. 

5.  ASSUMPTIONS  AND  LIMITATIONS:  None 

6.  ERROR  RETURNS:  None 

7.  LINKAGE  STATEMENT  AND  DESCRIPTION  OF  ARGUMENTS: 

CALL  MAXN  (X,  NX,  RSLT) 

entry  point  CALL  MINN  (X,  NX,  RLST) 

X  -  Array  frow  which  the  maximum  (minimum)  entry  is  to  be  found. 

NX  -  Integer,  number  of  entries  from  array  X  to  be  considered. 

RSLT  -  Maximum  (minimum)  entry  from  array  X. 

8.  PROGRAMS  CALLING  THIS  PROGRAM: 

MAXN:  S0RTR 

entry  point  MINN:  S0RTR 

9.  COMMON  VARIABLES  USED:  None 

10.  COMMON  VARIABLES  TO  BE  SET:  None 

11.  COMMON  VARIABLES  CHANGED:  None 

12.  PROGRAMS  USED  AND  DESCRIPTION  OF  LINKAGES:  None 

13.  NOTES  ON  METHODOLOGY:  If  entry  is  via  MAXN,  the  mode 

index  (M)  is  s'  -  to  2  indicating  the-"  the  maximum  value  of  the 
array  is  desired.  Entry  via  MINN  sets  M  to  1  indicating  that 
the  minimum  value  is  oesired  from  the  array.  At  statement 
1000  the  array  index  (MARK)  of  the  maximum  (or  minimum)  is 
initialized  to  one.  A  DO  loop  which  ranges  over  all  the  entries 
(NX)  considered  from  the  array  X  is  entered.  If  M  is  one,  statement 
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1030  is  used  to  determine  if  the  array  entry  being  tested 
(X(U)  is  less  than  the  currently  stored  minimum.  X(MARK).  tf 
X (L)  is  less  than  X(MARK),  control  transfers  to  statement 
10*11  where  the  array  index  (HARK)  for  the  minimum  array  entry 
is  set  to  the  present  entry  index  (L) .  If  M  .s  two,  statement 
1020  is  used  to  determine  ;f  the  array  entry  being 
tested  (X(L))  is  greater  than  the  currently  stored  maximum 
X(MARK).  If  X(L)  is  greater  than  X(MARK),  control  transfers 
to  statement  10*ll  where  the  array  index  (MARK)  for  the 
maximum  array  entry  is  set  to  the  present  entry  index  (L) . 

Upon  completion  of  the  above  loop  the  maximum  (M«2)  or 
minimum  (M*l)  value  for  the  array  is  stored  in  RSLT  and 
control  returned  to  the  calling  program. 
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1.  NAME:  MVC 

2.  TYPE  OF  PROGRAM:  SUBROUTINE,  LOGICAL 

3.  LANGUAGE:  FORTRAN  EXTENDED 

4.  PURPOSE:  Manipulate  character  strings  to  form  UTM  grid  square 
designators. 

5.  ASSUMPTIONS  AND  LIMITATIONS:  If  the  arguments  I0PTR  and  INPTR 
are  assigned  values  less  than  zero,  the  characters  will  not  be 
moved. 

6.  ERROR  RETURNS:  None 

7.  LINKAGE  STATEMENT  AND  DESCRIPTION  OF  ARGUMENTS: 

CALL  MVC  ( I CNT,  INPTR,  ZINADD,  I0PTR,  ZI0ADD) 

ICNT  -  Integer,  number  of  characters  to  be  moved  from  ZINADD 
to  ZI0ADD. 

INPTR  -  Integer,  index  of  the  first  character  to  be  moved  from 
ZINADD. 

ZINADD  -  Character,  2  word  array  (20  character)  from  which 
characters  will  be  moved. 

I0PTR  -  Integer,  index  of  the  location  in  ZI0ADD  where  the  first 
character  will  be  moved. 

ZI0ADD  -  Character,  2  word  array  (20  characters)  where  the 
characters  will  be  moved. 

8.  PROGRAMS  CALLING  THIS  PROGRAM:  ABSC0R,  CLN,  C00RT,  LABEL,  MAIN, 
MAP 

9.  COMMON  VARIABLES  USED:  None 

10.  COMMON  VARIABLES  TO  BE  SET:  None 

11.  COMMON  VARIABLES  CHANGED:  None 

12.  PROGRAMS  USED  AND  DESCRIPTION  OF  ARGUMENTS:  None 
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13.  NOTES  ON  METHODOLOGY:  Upon  entry  to  subroutine  MVC,  the  word 
location  (IW0RD)  in  the  array  ZI0ADD  where  the  characters  are 
to  be  moved  and  the  index  (ICHP0S)  of  the  first  character 
position  within  that  word  to  be  filled  are  calculated.  If 
ICHP0S  is  negative,  ICHP0S  is  set  to  9.  Next,  the  location 
(IW0D)  of  the  word  in  the  array  ZINADD  from  which  the  character 
string  is  to  be  moved  and  the  index  (ICHP0)  of  the  first 
character  position  to  be  moved  from  that  word  are  calculated. 

If  ICHP0  is  negative,  ICHP0  is  reset  to  9.  Next,  a  DO  loop 
is  entered  which  ranges  over  the  number  of  characters  to  be 
moved  (ICNT). 

Upon  entry  to  this  loop,  the  position  (ICHP0S)  in  word 
ZI0ADD  (W0RD)  where  the  first  character  will  be  moved  is  found 
by  adding  one  to  the  character  index  ICHP0S.  If  ICHP0S  lies 
within  word  ZI0ADD  (W0RD)  (ICHP0S  less  than  or  equal  10)  control 


transfers  to  statement  10.  If  ICHP0S  is  greater  than 
10,  the  character  position  is  reset  to  one  and  the  word 
location  IW0RD  incremented  by  one.  At  statement  10  the 
position  (ICH0)  of  the  first  character  to  be  moved  from  the 
word  ZINADD  (IW0R)  is  found  by  incrementing  ICH0  by  one. 

ICHP0  is  then  tested.  If  ICHP0  is  greater  than  10,  the  character 
position  lies  in  the  following  word  and  ICHP0  is  reset  to  one  and 
the  word  location  incremented  by  one.  At  statement  12,  TEMPI 
is  set  to  the  character  string  ZI0ADD  (IW0RD)  with  the  ICHP0S 
character  position  set  to  00  octal.  TEMP2  is  set  to  a  field 
of  blanks  with  the  ICHP0  character  position  set  to  the  ICHP0 
character  position  of  ZINADD  (IW0R).  TEMP2  is  then  shifted 
by  the  difference  in  position  of  the  two  characters  and  masked 
by  an  octal  7  in  character  position  ICHP0S.  This  result  is 
stored  in  TEMP2  ZI0ADD  (IW0RD)  is  then  the  0R  mask  of  TEMPI 
and  TEMP2.  The  loop  increments  ICHP0S  and  ICHP0  until  ICNT 
characters  have  been  moved  and  control  returns  to  the  calling 
program. 
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1.  NAME:  MyTIME 

2.  TYPE  OF  PROGRAM:  FUNCTION,  MATHEMATICAL 

3.  LANGUAGE  USED:  FORTRAN  EXTENDED 

k.  PURPOSE:  Calculate  accumulated  CP  time  in  units  of  lOOths  of  a 
second 

5.  ASSUMPTIONS  AND  LIMITATIONS:  None 

6.  ERROR  RETURNS:  None 

7.  LINKAGE  STATEMENT  AND  DESCRIPTION  OF  ARGUMENTS: 

I  «  MYTIME  (I) 

I  -  Integer,  accumulated  CP  time  (lOOths  seconds). 

8.  PROGRAMS  CALLING  THIS  PROGRAM:  JTIMER,  PP.ODMF 

9.  COMMON  VARIABLES  USED:  None 

10.  COMMON  VARIABLES  SET:  None 

11.  COMMON  VARIABLES  CHANGED:  None 

12.  PROGRAMS  USED  AND  DESCRIPTION  OF  LINKAGES:  None 

13.  NOTES  ON  METHODOLOGY:  Function  MYTIME  calculates  the  accumulated 

Control  Processor  (CP)  time  in  units  of  lOOths  of  a  second.  Upon 
entry  MYTIME  uses  the  function  SECOND  from  the  FORTRAN  EXTENDED 
library  to  find  the  accumulated  CP  time  in  second.  MYTIME  is 
calculated  by  multiplying  th'.s  dme  by  100  to  get  accumulated 
time  in  lOOths  of  a  second.  Control  is  then  returned  to  the 
calling  routine. 
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1.  NAME:  0NEPT 

2.  TYPE  OP  PROGRAM:  SUBROUTINE,  MATHEMATICAL 

3.  LANGUAGE:  FORTRAN  EXTENDED 

4.  PURPOSE:  Calculate  the  terrain  elevation  at  a  prespecified  grid 
point. 

5.  ASSUMPTIONS  AND  LIMITATIONS:  The  terrain  data  is  accessible  from 
the  direct  access  file  and  the  terrain  directory  has  been  previously 
stored. 

6.  ERROR  RETURNS:  A  RETURN  Al  will  be  executed  if  either  of  the  follow¬ 
ing  errors  are  detected: 

A.  Tiie  prespecified  point  coordinates  cannot  be  converted  to  UTM. 

B.  The  2  letter  UTM  grid  square  designator  for  the  square  of  interest 
cannot  be  converted  to  relative  ea.\ting  and  northing. 

C.  The  grid  square  containing  the  point  under  consideration  is  not 
on  the  direct  access  terrain  file. 

7.  LINKAGE  STATEMENT  AND  DESCRIPTION  OF  ARGUMENTS: 

CALL  0NEPT  (XS,  YS,  H) ,  RETURNS  (Al) 

XS  -  Relative  easting  of  the  grid  point  under  consideration. 

YS  -  Relative  northing  of  the  grid  point  under  consideration. 

H  -  MSL  altitude  of  the  terrain  at  the  grid  point  under  consideration. 

Al  -  Exit  which  is  taken  If  one  of  the  errors  described  in  6  (above) 
occurs . 

8.  PROGRAM  CALLING  THIS  PROGRAM:  ELEV 

9.  COMMON  VARIABLES  USED: 

A.  Z2BUF  -  I BUF 

B.  ZZFASC  -  L0DA,  LTDA,  N0DA 
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c.  ZZCV2  -  LTRS,  NBRS 
10.  COMMON  VARIABLES  TO  BE  SET: 

A.  ZZBUF  -  IBUF 

8.  ZZFASC  -  L0D.A,  LTDA,  N0DA 

C,  ZZCV2  -  LTRS 

H.  COMMON  VARIABLES  CHANGED: 

A.  ZZCV2  -  NBRS 

l2.  PROGRAMS  USED  AND  DESCRIPTION  OF  LINKAGES: 

A.  SUBROUTINES 

(2)  CALL  cm  (X.  V,  D.  RETURNS  (Ai)  (entry  point  to  C00RT)- 
Converts  the  re.ative  eesting  end  northing  of  the  Ro¬ 
dder  considerotion  to  UTH  coordinates  (2  letters, 
x  .  (tetative  easting  of  the  point  under  consideration, 
y  -  Relative  northing  of  the  point  under  consideration. 

L  -  gunny  variable  not  used  for  this  entry  point. 

A,  -  Control  is  transferred  to  this  statement  number  if 
an  error  is  detected  in  C1TA. 

(31  CALL  CITN  (X,  Y,  L) ,  RETURNS  (Al)  (entry  point  to  C00RT)  - 
Converts  the  UTM  coordinates  of  a  point  to  re  at.ve 

easting  and  northing. 

X  -  Relative  easting  of  the  point  under  consideration, 
y  -  Relative  northing  of  the  point  under  consideration. 

L  -  Dumny  variable  not  used. 

Al  -  Control  is  transferred  to  this  statement  number  if 
an  error  is  detected  in  ClTN. 
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13.  NOTES  ON  METHODOLOGY:  Upon  entry  to  0NEPT ,  the  lateral  data  rate, 
RATE,  is  set  to  500  meters.  C1TA  is  then  called  to  calculate  the  UTM 
coordinates  of  the  location  (XS,Y$).  A  DO  loop  is  then  used 
to  determine  if  the  grid  square  (LTRS)  which  contains  point  (X5,  YS) 
is  available  from  the  direct  access  terrain  file,  if  this  grid 
square  is  not  on  the  terrain  file,  control  transfers  to  statement 
9000  where  a  RETURN  Al  ;s  executed.  If  the  grid  square  is 
available  from  the  direct  access  file,  control  transfers  to  statement 
1300  where  the  address  of  the  first  'e-rain  record  for  this  square 
is  stored  in  local  variable  L0C.  C1TN  is  then  called  to  calculate 
the  relative  easting  and  northing  of  the  southwest  corner  of  this 
grid  square. 

A  DO  loop  is  then  entered.  This  loop  ranges  over  the  24  records 
containing  the  terrain  elevations  for  this  100  km  grid  square.  Each 
record  covers  a  rectangular  area  where  the  sides  are  defined  by 
the  boundaries  XI,  (western),  X2  (eastern),  Y1  (southern)  and 
Y2  (northern).  A  test  is  then  performed  to  determine  if  the 
point  (XS,  YS)  is  in  this  rectangle.  If  the  point  is  not  within 
this  record  area,  control  transfers  to  the  end  of  the  loop.  The 
record  number  is  incremented  by  one  and  the  next  record  area  is 
tested.  When  the  record  which  covers  (XS,  YS)  has  been  found, 
then  the  point  number  within  the  record  and  the  record  address 
IR  is  not  equai  to  LASL0C  (1),  the  required  record  is  read 
from  the  direct  access  file  and  stored  in  IBUF.  The  MSL 
altitude  is  then  calculated  by  multiplying  the  altitude  IBUF  (IB) 
by  the  vertical  data  rate  (10  meters).  Control  then  returns  to 
the  calling  program. 
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1 .  NAME :  PACK 

2.  TYPE  OF  PROGRAM:  SUBROUTINE,  OUTPUT 

3.  LANGUAGE:  FORTRAN  EXTENDED 

k.  PURPOSE:  Punch  a  track  Identifier,  track  point  identifier 
and  logical  flag  on  file  9. 

5.  ASSUMPTIONS  AND  LIMITATIONS:  None 

6.  ERROR  RETURNS:  Hone 

7.  LINKAGE  STATEMENT  AND  DESCRIPTION  OF  ARGUMENTS: 

Call  PACK  (I ,  ISITE,  ZWENG) 

I  -  Integer,  track  index. 

ISITE  -  Integer,  site  index. 

ZWENG  -  Logical,  site-track  engageability  flag. 

8.  PROGRAMS  CALLING  THIS  PROGRAM:  MAIN 

9.  COMMON  VARIABLES  USED:  None 

10.  COMMON  VARIABLES  TO  BE  SET:  None 

11.  COMMON  VARIABLES  CHANGED:  None 

12.  PROGRAMS  USED  AND  DESCRIPTION  OF  LINKAGES:  None 

13-  NOTES  ON  METHODOLOGY:  Each  time  PACK  is  called,  a  record 

is  punched  on  file  9  containing  the  track  identifier,  the 
track  point  identifier  and  the  engagement  flag.  Control 
is  then  returned  to  the  calling  program. 
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1.  NAME:  PMAP 

2.  TYPE  OF  PROGRAM:  SUBROUTINE,  LOGICAL 

3.  LANGUAGE:  FORTRAN  EXTENDED 

4.  PURPOSE:  Control  the  input  of  the  terrain  data  from  tape  and 
initialize  program  execution. 

5.  ASSUMPTIONS  AND  LIMITATIONS:  File  7  contains  the  terrain  data 
tape. 

6.  ERROR  RETURNS:  None 

7.  LINKAGE  STATEMENT  AND  DESCRIPTION  OF  ARGUMENTS:  The  linkage  to 
PMAP  is  from  the  SC0PE  operating  system. 

8.  PROGRAMS  CALLING  THIS  PROGRAM:  None 

9.  COMMON  VARIABLES  USED: 

A.  2ZRAI  -  IMAX 

B.  ZZTRKD  -  NDWfRD 

10.  COMMON  VARIABLES  TO  BE  SET:  None 

11.  COMMON  VARIABLES  CHANGED: 

A.  ZZRAI  -  IMAX 

B.  ZZTRUD  -  NDW0RD 

12.  PROGRAMS  USED  AND  DESCRIPTION  OF  LINKAGES: 

A.  SUBROUTINES 

(I)  CALL  MAIN,  RETURNS  (Al)  -  Controls  input  of  data  cards, 
checking  of  input  data  and  plotting  of  maps. 

Al  -  Control  is  transferred  to  this  statement  number  if 
an  error  is  detected  ir.  HAIN. 
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(2)  CALL  RAW  (TIN,  T0UT) ,  RETURNS  (Al)  -  Read  terrain  data 
from  a  magnetic  tape  (file  TIN)  and  write  this  data  on 
the  direct  access  file  T0UT. 

TIN  -  Integer,  file  number  of  tape  to  be  read. 

T0UT  -  Integer,  file  number  of  direct  access  device  where 
the  terrain  data  will  be  stored. 

Ai  -  Control  is  transferred  to  this  statement  number  if 
an  error  is  detected  in  RAW. 

13.  NOTES  ON  METHODOLOGY:  Upon  entry  to  PMPP,  NDW0RD  is  set  to  3000, 

KINN  is  set  to  7  and  K0UT  is  set  to  4.  Si" IP  and  NSQ.UAR  are  then 
read  from  the  terrain  tape  (file  7).  The  maximum  number  of  terrain 
records  .(IMAX)  to  be  read  from  the  terrain  file  is  then  calculated. 

RAW  is  then  called  to  read  the  terrain  data  from  the  tape  (file  7) 
and  store  this  terrain  data  on  the  direct  access  device  (file  4). 

MAIN  is  "hen  called  to  provide  the  logical  sequencing  required  for 
successful  PMAP  execution.  Control  then  returns  to  the  SC0PE  operating 
system. 
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1.  NAME:  RAW 

2.  TYPE  OF  PROGRAM:  SUBROUTINE,  INPUT/jfUTPUT 

3.  LANGUAGE:  FORTRAN  EXTENDED 

4.  PURPOSE:  Read  the  terrain  data  file  from  type,  create  a  directory, 
and  write  the  directory  and  terrain  data  on  a  direct  access  storage 
file. 

5.  ASSUMPTIONS  AND  LIMITATIONS:  The  tape  (file  10)  has  the  directory  and 
terrain  data  In  the  proper  sequence.  The  direct  access  storage  file 
(file  11)  must  be  made  available. 

6.  ERROR  RETURNS:  A  RETURN  Al  will  occur  if  the  following  error  Is 
detected: 

A.  The  number  of  records  read  exceeds  the  maximum  I MAX. 

7.  LINKAGE  STATEMENTS  AND  DESCRIPTION  OF  ARGUMENTS: 

CALL  RAW  (TIN,  TOUT),  RETURN  (Al) 

TIN  -  Number  of  the  input  tape  file  (10). 

TOUT  -  Number  of  the  output  direct  access  file  (11). 

Al  -  Exit  taken  If  the  condition  listed  under  6  (above)  occurs. 

8.  PROGRAMS  CALLING  THIS  PROGRAM:  PMAP 

9.  COMMON  VARIABLES  USED: 

A.  2ZBUF  -  IBUF 

B.  ZZRAI  -  I  MAX,  IRAN 

10.  COMMON  VARIABLES  TO  BE  SET: 

A.  ZZRAI  -  I  MAX 

11.  COMMON  VARIABLES  CHANGED: 

A.  ZZBUF  -  IBUF 

B.  ZZRAI  -  IRAN 
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12.  PROGRAHS  USED  AND  DESCRIPTION  OP  LINKAGES:  None. 

13.  NOTES  ON  METHODOLOGY:  RAW  reads  the  terrain  data  from  tape  (file 
TIN)  and  writes  these  data  on  the  direct  access  device  (file  T0UT). 
Upon  entry  to  RAW  a  REWIND  is  executed  for  both  files  to  provide  the 
proper  starting  reference  for  these  files.  Sixty-two  words  containing 
the  grid  zone  identifier,  grid  square  identifiers  and  terrain  data 
locations  are  read  from  tape  (file  TIN),  stored  in  IBUF,  and  then 
written  on  the  direct  access  device  (file  T0UT). 

Next  the  record  counter  (N0  REC)  is  set  to  one  and  a  loop 
entered.  This  loop  increments  the  record  counter,  reads  a  record 
(1717  or  1515  terrain  data  points)  from  tape  and  writes  this  data  on 
the  direct  access  device  (file  T0UT).  This  loop  reads  records  from 
the  tape  and  writes  them  on  the  direct  access  device  until  an  end  of 
file  mark  is  encountered  or  the  maximum  number  of  records  (IMAX)  is 
exceeded.  If  IMAX  is  exceeded,  control  is  transferred  to  statement 
AO  where  an  error  message  is  printed,  a  REWIND  TIN  is  executed,  and 
an  error  exit  (RETURN  Al)  taken.  If  an  end  of  file  is  encountered, 
control  transfers  to  statement  80  where  a  REWIND  TIN  Is  executed  and 
a  normal  return  taken. 
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1.  NAME:  SEARCH 

2.  TYPE  OF  PROGRAM:  SUBROUTINE,  LOGICAL 

3.  LANGUAGE:  FORTRAN  EXTENDED 

4.  PURPOSE:  Find  the  targeted  track  point  for  a  specified  site. 

5.  ASSUMPTIONS  AND  LIMITATIONS:  None 

6.  ERROR  RETURNS:  None 

7.  LINKAGE  STATEMENT  AND  DESCRIPTION  OF  ARGUMENTS: 

CALL  SEARCH  (I SITE) 

ISITE  -  Integer,  specified  site  index. 

8.  PROGRAMS  CALLING  THIS  PROGRAM:  MAIN 

9.  COMMON  VARIABLES  USED:  None 

10.  COMMON  VARIABLES  TO  BE  SET:  None 

11.  COMMON  VARIABLES  CHANGED:  None 

12.  PROGRAMS  USED  AND  DESCRIPTION  OF  LINKAGES:  Hone 

13.  NOTES  ON  METHODOLOGY:  Upon  entry  to  subroutine  SEARCH,  a  REWIND 

is  executed  on  file  9  to  assure  that  this  file  starts  at  the  proper 
position.  The  first  record  is  read  from  file  9  at  statement  450. 

If  an  end  file  is  encountered,  control  transfers  to  statement 
999  where  a  normal  return  is  taken  to  the  calling  program.  At 
statement  970,  a  check  is  performed  to  determine  if  an  1/0 
parity  error  was  detected  during  this  read.  If  an  1/0 
parity  error  was  detected,  control  transfers  to  statement  990  where 
an  error  message  is  printed  and  control  returns  to  the  calling 
program.  If  an  end  of  file  is  not  encountered  and  an  1/0  parity 
error  is  not  detected,  the  track  point  identifier  (I DTP)  read  from 
this  record  is  tested  against  the  site  identifier  (ISITE)  specified 
in  the  argument  list.  If  these  identifiers  are  not  egual,  control 
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transfers  to  statement  450  where  the  next  record  will  be  read. 

If  these  identifiers  are  equal,  the  engagement  flag  ZWENC  is 
tested.  If  the  engagement  flag  is  TRUE,  a  targeted-by-track 
message  is  printed  and  control  returns  to  the  catling  subroutine. 
Otherwise,  control  returns  to  statement  450  where  the  next 
record  will  be  read. 
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1.  NAME:  SETBPT 

2.  TYPE  OF  PROGRAM:  SUBROUTINE,  LOGICAL 
.  LANGUAGE:  FORTRAN  EXTENDED 

4.  PURPOSE:  Calculate  the  block  points  address  for  a  specified 
track  and  point. 

5.  ASSUMPTIONS  AND  LIMITATIONS:  None 

6.  ERROR  RETURNS:  None 

7.  LINKAGE  STATEMENT  AND  DESCRIPTION  OF  ARGUMENTS: 

CALL  SilfBPT  (NTK,  I,  IW0RD) 

NTK  -  Index  for  the  desired  track. 

I  -  Index  of  the  desired  end  point. 

IW0RD  -  Block  pointer  address  for  desired  point. 

8.  PROGRAMS  CALLING  THIS  PROGRAM:  JTRK,  TTRK,  ZCKPTI 

9.  COMMON  VARIABLES  USED: 

A.  ZZTRKD  *  I  DATA  (equi valenced  to  TKDATA) 

10.  COMMON  VARIABLES  TO  BE  SET: 

A.  ZZTRKD  -  I  DATA  (equi valenced  to  TKDATA) 

11.  COMMON  VARIABLES  CHANGED:  None 

12.  PROGRAMS  USED  AND  DESCRIPTION  OF  LINKAGES:  None 

13.  NOTES  ON  METHODOLOGY:  Upon  entry  to  SETBPT  the  location  of  the 
word  ( iWL0c)  in  the  track  data  array  containing  this  block  pointer 
address  is  calculated  for  this  track.  The  block  number  (NBLK)  is 
then  determined  via  a  mask  expression  and  shift  of  the  word  found  at 
location  1WL0C.  The  block  pointer  word  ( IW0RD)  is  then  determined  from 
NBLK  and  the  point  index  I.  If  the  debug  flag  is  set  to  TRUE, 

the  track  data  array  is  dumped.  Control  then  returns  to  the  calling 
program. 
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1 .  NAME:  S0RTR 

2.  TYPE  OF  PROGRAM:  SUBROUTINE,  LOGICAL 

3.  LANGUAGE:  FORTRAN  EXTENDED 

k.  PURPOSE:  Perform  either  ascending  or  descending  sort  on  an  array 
and  output  maximums  and  minimums. 

5.  ASSUMPTIONS  AND  LIMITATIONS:  None 

6.  ERROR  RETURNS:  None 

7.  LINKAGE  STATEMENT  AND  DESCRIPTION  OF  ARGUMENTS: 

CALL  S0RTR  (X,Y,NX,MAXNX,NY,M0DE) 

X  -  Array  for  which  an  ascending  sort,  descending  sort  or  maximum 
and  minimum  of  the  array  entries  will  be  calculated. 

Y  -  Array  whose  entries  will  be  moved  in  an  order  corresponding 
to  the  sort  performed  on  X. 

NX  -  Integer,  number  of  entries  from  the  array  X  to  be  available 
for  sort. 

MAXNX  -  Integer,  maximum  number  of  entries  upon  which  a  sort  is 
to  be  performed. 

NY  -  Integer,  numcer  of  entries  in  array  Y  which  will  be  moved 
for  each  corresponding  move  in  the  X  array. 

MODE  -  Integer,  mode  of  operation  for  this  routine.  1  implies 
descend : ng  sort j  2  implies  ascending  sort |  3  impli**^ 
maximum  and  minimum  entries  only  are  desired. 

8.  PROGRAM  CALLING  THIS  PROGRAM:  MAIN 
9-  COMMON  VARIABLES  USED: 

A.  ZZMAPP  -  XMAX,  YMAX,  XMIN,  YMIN 
10.  COMMON  VARIABLES  TO  BE  SET:  None 
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U.  COMMON  VARIABLES  CHANGED: 

A.  ZZMAPP  “  XMAX,  YMAX,  XMIN,  YMIN 
12.  PROGRAMS  USED  AND  DESCRIPTION  OF  LINKAGES: 

A.  SUBROUTINES 

(I)  CALL  MAXN(X,NX,RSLT)  -  Find  the  maxima  ontry  In  array  X. 

X  -  Array  from  which  the  maximum  entry  Is  to  be  found. 

NX  -  Integer,  number  of  entries  from  array  X  to  be 
considered. 

RSLT  -  Maximum  entry  from  array  X. 

(2)  CALL  M1NM(X,NX,RSLT)  (entry  point  to  MAXN)  -  Find  the 
the  minimum  entry  In  array  X. 

X  -  Array  from  which  the  maximum  entry  is  to  be  found. 

NX  -  Integer,  number  of  entries  from  array  X  to  be 
considered. 

RSLT  -  Maximum  entry  from  array  X. 

13.  NOTES  ON  METHODOLOGY:  Upon  entry  to  S0RTR ,  h  «s  set  to  the 

number  of  array  entries  to  be  sorted  (NX)  minus  one-  At  statement 
1010  the  sort  DO  loop  is  entered.  First  the  lower  (Ll)  and  upper 
(LL)  indices  which  determine  the  entries  which  require  sort  are 
determined.  The  lower  index  is  set  to  the  DO  loop  Index  K  and 
the  upper  Index  is  set  to  NX.  The  array  index  MARK  is  set  to  the  00 
loop  index  K.  At  statement  1012,  a  DO  loop  which  ranges  over  the 
entries  (Ll  to  LL)  to  be  sorted  from  the  array  X  is  entered.  If 
MiOE  is  one,  statement  1020  is  used  to  determine  if  the  array 
entry  being  tested  (X(D)  is  greater  than  the  currently  stores 
maximum,  X(MARK) ,  If  X(L)  is  greater  than  X(MARK),  control  transfers 
to  statement  1041  where  the  array  Index  (MARK)  for  the 
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maximum  array  entry  is  set  to  the  present  entry  index  (L) .  If 
M0DE  is  two,  statement  1030  is  used  to  determine  if  the 
array  entry  being  tested  (X(U)  is  less  than  the  currently  stored 
minimum,  X (MARK) .  If  X(L)  is  less  than  X(MARK),  control  transfers 
to  statement.  1041  where  the  array  index  (MARK)  for  the 
minimum  array  entry  is  set  to  the  present  entry  index  (L) .  If 
M0DE  is  3,  no  sorts  are  required  and  control  transfers  to  statement 
1052.  Following  statement  1049  the  maximum  (M0DE*1)  or 
minimum  (M0D£=2)  array  entry  is  stored  in  the  array  location  K 
and  the  value  of  the  array  previously  stored  in  array  location  K 
is  stored  in  location  MARK.  This  procedure  stores  the  values 
X(MARK)  in  descending  order  (M0DE=1)  or  ascending  order  (M0DE=2) 
as  required. 

If  NY  is  greater  than  zero,  NY  entries  from  the  array  Y  wi 1 1 
be  switched  corresponding  to  the  switch  performed  in  the  X  array. 

When  the  sort  is  complete  or  IM0DE  is  3,  control  transfers 
to  statement  1052.  Following  statement  1052  the  maximum  (XMAX)  and 
minimum  (XMIN)  entries  in  array  X  are  found.  If  NY  is  less  than 
or  equal  to  zero,  control  returns  to  the  calling  program. 

Otherwise  the  maximum  (YMAX)  and  minimum  (YMIN)  entries  for  the 
array  Y  are  calculated.  Control  then  returns  to  the  calling 
program. 
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1.  NAME:  STPNT  (other  entry  points:  STRK) 

2.  TYPE  OF  PROGRAM:  SUBROUTINE,  STORAGE 

3.  LANGUAGE:  FORTRAN  EXTENDED 

4.  PURPOSE:  Store  track  data  in  the  track  data  table. 

5.  ASSUMPTIONS  AND  LIMITATIONS:  STRK  must  be  called  to  calculate 
track  point  initial  address  prior  to  point  storage. 

6.  ERROR  RETURNS:  A  error  message  will  be  printed  and  a  RETURN  El 
executed  if  track  data  table  overflows. 

7.  LINKAGE  STATEMENT  AND  DESCRIPTION  OF  LINKAGES: 

CALL  STPNT  (X,  Y,  Z,  T,  J),  RETURNS  (El) 

entry  point  CALL  STRK  (DUM,  DUM,  DUM,  DUM.NDUM),  RETURNS  (El) 

X  -  Relative  easting  of  the  track  point  under  consideration. 

Y  -  Relative  northing  of  the  track  point  under  consideration. 

Z  -  Altitude  of  the  track  point  under  consideration. 

T  -  Time  required  to  reach  the  point  under  consideration. 

J  -  Check  point  code  for  the  point  under  consideration. 

DUM  -  Dummy  variable  not  used. 

NDUM  -  Dummy  variable  not  used. 

El  -  Exit  taken  if  the  condition  listed  under  6  (above)  occurs. 

8.  PROGRAMS  CALLING  THIS  PROGRAM: 
entry  point  STPNT:  MAIN 
entry  point  STRK:  MAIN 

9.  COMMON  VARIABLES  USED: 

A.  ZZTRKD  -  NDW0RD,  TKDATA  (equi valenced  to  I  DATA) 

10.  COMMON  VARIABLES  TO  BE  SET: 

A.  ZZTRKD  -  NDW{*R0 

11.  COMMON  VARIABLES  CHANGED: 

A.  ZZTRKD  -  TKDATA  (equi  %'alenced  to  I  DATA) 
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12.  PROGRAMS  USED  AND  DESCRIPTION  OF  LINKAGES:  None 

13.  NOTES  ON  METHODOLOGY:  Upon  entry  to  STPNT,  the  path  leg  index 

(LEG)  is  incremented  and  the  address  of  the  first  word  of  storage 

for  this  leg  is  calculated.  If  this  storage  area  exceeds  the  max¬ 

imum  word  address  (NDWORD) ,  control  transfers  to  statement  9000 
where  an  overflow  message  is  printed  and  control  transfers  to 

the  calling  program  via  error  RETURN  El.  If  table  overflow  does 
not  occur,  the  time  to  reach  this  leg  endpoint  is  stored  in  the 

track  data  array  (TKDATA)  at  address  IWORDA.  The  easting  of  this 

point  is  then  stored  in  the  first  40  bits  of  the  next  word  and  the 
and  the  corresponding  track  altitude  is  stored  in  the  last  20 
bits  of  this  word.  The  address  is  then  incremented  and  the 
points  northing  stored  in  the  first  bits  of  the  .iext  word  and  the 
checkpoint  code  (MC)  stored  in  the  last  20  bits  of  this  word. 

Control  then  transfers  to  statement  5000  where  the  track  array 
will  be  dumped  if  the  DEBUG  flag  ( DB)  is  TRUE.  Control  is  then 
returned  to  the  calling  program. 

Upon  entry  via  STRK,  the  total  leg  counter  (LEGSUM)  is  incre¬ 
mented  by  the  number  of  the  legs  processed  for  the  previous  track. 
The  location  counter  (NL0C)  is  then  incremented  by  one  and  the 
address  (IFW0RD)  of  the  first  word  in  this  block  is  calculated. 

Then  the  block  counter  (NBLK)  is  set.  If  this  is  the  first  call 
to  STRK,  the  track  data  array  is  initialized  to  zero.  At  state¬ 
ment  A000,  the  location  of  the  block  data  pointer  address  is 
calculated.  The  block  data  address  is  then  stored  in  the  data 
array.  Control  then  returns  to  the  calling  program. 
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J.  NAME:  THYME  (other  entry  points:  TZT) 

2.  TYPE  OF  PROGRAM:  SUBROUTINE,  MATHEMATICAL 

3.  LANGUAGE:  FORTRAN  EXTENDED 

k.  PURPOSE:  Calculate  and  print  the  program  accumulated  time  and  incre¬ 
mental  time  since  the  last  time  printout  via  THYME. 

3.  ASSUMPTIONS  AND  LIMITATIONS:  None 

6.  ERROR  RETURNS:  None 

7.  LINKAGE  STATEMENT  AND  DESCRIPTION  OF  ARGUMENTS: 

CALL  THYME  no  arguments 

entry  point  CALL  TZT  no  arguments 

8.  PROGRAMS  CALLING  THIS  PROGRAM: 
entry  point  THYME: 

entry  point  TZT:  MAIN,  OUPINP,  PR0DMF,  PSIT,  PSYS 

9.  COMMON  VARIABLES  USED:  None 

10.  COMMON  VARIABLES  SET:  None 

11.  COMMON  VARIABLES  CHANGED:  None 

12.  PROGRAMS  USED  AND  DESCRIPTION  OF  LINKAGES: 

A.  FUNCTIONS 

(I)  TIME  =  MYTIME(M).  HYTIME  is  used  to  calculate  the  accumulated 
CP  time  in  units  of  lOOths  of  a  second. 

M  -  Integer,  dummy  variable  not  used. 

13.  NOTES  ON  METHODOLOGY:  Subroutine  THYME  is  used  to  calculate  accu¬ 
mulated  CP  time  and  incremental  CP  times  and  print  a  two  line  time 
report.  Upon  entry  to  THYME  a  printer  line  is  skipped.  Entry  via 
ENTRY  TZT  first  calculates  TIME,  the  accumulated  CP  time  in  seconds, 
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using  function  MYTIME.  A  printout  of  the  accumulated  time,  TIME, 
is  then  provided.  If  ITIM,  the  accumulated  time  (stored  from  pre¬ 
vious  call),  is  less  than  .001  seconds  control  transfers  to  state¬ 
ment  1050  where  TIMI  is  set  to  the  current  accumulated  time  and  a 
normal  exit  to  the  calling  subroutine  taken.  If  ITIM  is  greater 
than  or  equal  to  .001,  the  incremental  CP  time  is  calculated  from 
the  difference  between  the  current  time,  TIME,  and  previous  time, 
ITIM.  A  printout  of  this  incremental  time  is  provided,  ITIM  is  set 
to  TIME  and  control  is  returned  to  the  calling  program. 
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1.  NAME:  TT0RNG 

2.  TYPE  OF  PROGRAM:  SUBROUTINE,  MATHEMATICAL 

3.  LANGUAGE:  FORTRAN  EXTENDED 

4.  PURPOSE:  Calculate  the  time  a  specified  track  enters  a  specified 
site  fire  volume. 

5.  ASSUMPTIONS  AND  LIMITATIONS:  A  entry  is  recorded  only  if  this  time 
is  after  TI2Z. 

6.  ERROR  RETURNS:  None 

7.  LINKAGE  STATEMENT  AND  DESCRIPTION  OF  ARGUMENTS: 

CALL  TT0RNG  (RAG IN,  TIZZ,  TIME,  NF,  HTK) 

RAGIN  -  Range  of  interest  for  the  site  under  consideration. 

TIZZ  -  Earliest  time  to  be  considered  for  entry. 

TIME  -  Time  a  specified  track  enters  a  specified  site  fire  volume. 
NF  -  Integer,  index  of  the  site  under  consideration. 

NTK  -  Integer,  index  of  the  track  under  consideration. 

8.  PROGRAMS  CALLING  THIS  PROGRAM:  MAIN 

9.  COMMON  VARIABLES  USED: 

A.  ZZTRK1  -  MNP 

B.  ZZSITE  -  XFU,  YFU,  ZFU 

10.  COMMON  VARIABLES  TO  BE  SET: 

A.  ZZTRKI  -  MNP 

B.  ZZSITE  -  XFU,  YFU,  ZFU 

11.  COMMON  VARIABLES  CHANGED:  None 

12.  PROGRAMS  USED  AND  DESCRIPTION  OF  LINKAGES: 

A.  FUNCTIONS 

(I)  X  =  TTRK  (NTK,l)-  Finds  the  time  for  the  specified  track 
and  point  indices. 
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NIK  -  Index  for  the  specified  track. 

I  -  Index  for  the  specified  point. 

(2)  X.*  XTRK  (NTK,  l)(entry  point  to  TTRK)  -  Finds  the  easting 

for  the  specified  track  and  point  indices. 

NTK  -  Index  for  the  specified  track. 


I  -  Index  for  the  specified  point. 

(3)  X_  *  YTRK  (NTR,  I)  (entry  point  tc  TTP,K)-Finds  the  northing 
for  the  specified  track  and  point  indices. 


NTK  -  Index  for  the  specified  track. 

I  -  Index  for  the  specified  point. 

(4)  X^ *  ZTRK  (NTK,  l)  (entry  point  to  TTRK)  -  Finds  the  altitude 
for  the  focified  'rack  and  point  indices. 


NTK  -  Index  for  the  specified  track. 
I  -  Index  for  the  specified  point. 


13.  NOTES  ON  METHODOLOGY:  UDon  entry  to  TT<RNG,  the  site  location  is 
stored  in  local  variables  XFUT,  YFUT,  AND  ZFUT.  The  square  of 
the  site  range  of  interest  is  calculated,  the  number  of  points 
on  the  track  is  stored  locally,  and  the  point  counter  ( 1 )  for  this 
track  is  initialized  to  1. 

At  statement  1005,  a  loop  is  entered  which  finds  the  time 
of  intersection  for  a  specified  track  with  a  specified  site 
fire  volume.  First,  the  point  counter  is  indexed  and  a  test 
performed  to  see  if  the  maximum  number  of  points  for  this 
track  has  been  exceeded.  If  the  maximum  number  of  points 
have  been  exceeded  (l>MN),  control  transfers  to  statement 
1036  where  the  time  of  intersection  Is  set  equal  to  -l  and 
control  returns  to  the  calling  program.  If  the  number  of 
points  has  not  been  exceeded,  a  test  is  performed  to  see  if 
this  leg  end  point  is  reached  before  the  specified  initial 
time  (TIZZ).  If  this  leg  end  point  is  reached  by  the 
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vehicle  before  the  specified  initial  time,  control  transfers  to 
5  statement  1005  whee  the  point  index  is  incremented  in 

t 

1  preparation  for  trie  next  leg  to  be  tested.  Next,  the  square  of 

the  range  from  the  site  location  to  the  leg  starting  point  and 
i  end  point  (ANS,CNS)  are  calculated.  Tire  is  set  equal  to  the 

leg  initialization  time  and  the  range  to  the  starting  point 
is  tested.  If  this  range  to  the  leg  starting  point  is  less 
than  the  range  of  interest,  control  returns  to  the  calling  pro- 
;  cram.  If  this  point  lie-,  'utside  the  range  of  interest,  the 

leg  extension  variable  DISCRM  is  calculated.  If  DISCRM  is  less 
than  zero,  this  leg  does  not  intersect  the  fire  volume  and  control 
transfers  to  statement  1005.  If  this  leg  intersects  the  fire 
volume,  the  entry  time  and  exit  time  for  this  leg  (or  leg  extension) 
are  calculated.  If  this  time  satisfies  the  desired  constraints, 
the  entry  time  is  stored  in  TIKE  rnd  control  then  transfers  to  the 
calling  program. 
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1.  NAME:  TTRK  (other  entry  points:  XTRK,  YTRK,  7.TRK) 

2.  TYPE  OF  PROGRAM:  FUNCTION,  LOGICAL 

3.  LANGUAGE:  FORTRAN  EXTENDED 

4.  PURPOSE:  Find  values  of  time,  easting,  northing,  and  altitude 
from  the  track  array. 

5.  ASSUMPTIONS  AND  LIMITATIONS:  None 

6.  ERROR  RETURNS:  None 

7.  LINKAGE  STATEMENT  AND  DESCRIPTION  OF  ARGUMENTS: 
entry  point  £  =  TTRK(NTK.l) 

entry  point  £  =  XTRK(NTK.l) 
entry  point  £  =  YTRK(NTK.l) 
entry  point  £=  ZTRK(NTK,l) 

NTK  -  Integer,  index  of  the  path  of  interest. 

I  -  integer,  index  of  the  point  of  interest. 

8.  PROGRAMS  CALLING  THIS  PROGRAM: 

TTRK:  MAIN,  TT0RNG 

entry  point  XTRK:  MAIN,  TT0RNG 
entry  point  YTRK:  MAIN  TT0RNG 
entry  point  ZTRK:  MAIN,  TT0RNG 

9.  COMMON  VARIABLES  USED: 

A.  ZZTRK1  -  MNP 

B.  ZZSITE  -  aFU,  YFU,  ZFU 

10.  COMMON  VARIABLES  TO  BE  SET: 

A.  ZZTRK1  -  MNP 

B.  ZZSITE  -  ZFU,  YFU,  ZFU 

11.  COMMON  VARIABLES  CHANGED:  None 
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12.  PROGRAMS  USED  AND  DESCRIPTION  OF  LINKAGES: 

A.  SUBROUTINES 

CALL  SETBPT  (NTK,  I , IW0RD; 

NTK  -  Index  of  the  desired  track. 

I  -  Index  of  the  desired  point. 

IW0RD  -  Block  pointer  address  for  the  desired  point. 

13.  NOTES  ON  METHODOLOGY:  Function  TTRKwith  entry  points  XTRK, 

YTRK,  end  ZTRK,  is  used  to  find  the  time,  easting  and  northing 
for  a  specified  track  point.  Upon  entry  via  TTRK,  the  block 

pointer  word  (IW0RD)  is  calculated  hv  a  nail  to  subroutine  SETBPT. 
The  time  is  then  stored  in  variable  TTRK  and  control  transfers 
to  statement  6000.  if  entry  is  via  entry  point  XTRK, 
the  address  index  (JA)  is  initialized  to  1  and  control  transfers 
to  statement  A000.  If  entry  is  via  entry  point  YTRK,  the 
addreu  index  is  initialized  to  2  and  control  transfers  to 
statement  A000.  At  statement  1*000,  SETBPT  is  called  to 
calculate  the  block  pointer  address  for  this  track  point.  The 
easting  (JA=l)  or  northing  (JA=2)  is  then  stored  and  control  trans 
fers  to  statement  6000.  If  entry  is  via  entry  point  ZTRK, 

SETBPT  is  called  to  find  the  block  pointer  address  for  this 
point.  The  altitude  at  this  point  is  then  calculated  at  state¬ 
ment  5000.  The  debug  flag  ( DB)  is  then  tested.  If  the 
debug  flag  is  TRUE,  PDUMP  is  called  to  dump  the  track  data  array 
and  control  returns  to  the  calling  program. 
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NAME:  ZCKPT1  (other  entry  points:  ZCKPT2,  ZCKPT3,  ZCKPT4) 

TYPE  OF  PROGRAM:  FUNCTION,  LOGICAL 
LANGUAGE:  FORTRAN  EXTENDED 

PURPOSE:  Set  the  logical  check  point  flap  correspondinq  to  the 
appropriate  bi t. 

ASSUMPTIONS  AND  LIMITATIONS:  None 
ERROR  RETURNS:  None 

LINKAGE  STATEMENT  AND  DESCRIPTION  OF  ARGUMENTS: 

entry  point  «*  ZCKPTI (NTK, I ) 
entry  point  X_  =  ZCKPT2 (NTK,  I ) 
entry  point  YL  =  ZCKPT3 (NTK,  I ) 
entry  point  X_  =  ZCKPTA{NTK, I ) 

All  use  the  same  arguments. 

NTK  -  Integer,  Index  of  the  path  of  interest. 

I  -  Integer,  Index  of  the  point  of  interest. 

PROGRAMS  CALLING  THIS  PROGRAM:  None 

COMMON  VARIABLES  USED: 

A.  ZZTRKD  -  TKDATA  (equivalenced  to  I  DATA) 

COMMON  VARIABLES  TO  BE  SET: 

A.  ZZTRKD  -  TKDATA  (equivalenced  to  I  DATA) 

COMMON  VARIABLES  CHANGED:  None 

PROGRAMS  USED  AND  DESCRIPTION  OF  LINKAGES: 

A.  SUBROUTINES: 

CALL  SETBPT  (NTK,  I,  IWjSRD) 

NTK  -  Index  of  the  desired  track. 

I  ■  Index  of  the  desired  point. 

IW0RD  -  Block  pointer  address  for  the  desired  point. 

NOTES  ON  METHODOLOGY:  Entry  to  this  function  via  entry  points 
ZCKPTI ,  ZCKPT2,  ZCKPT3,  and  ZCKPTA  specifies  the  checkpoint  bit 
index  for  this  track  point.  Upon  ertry  the  first,  second,  third, 
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or  fourth  bit  of  the  check  point  code  will  be  selected  for  test. 
A*:  statement  5000,  the  block  pointer  address  for  this 
track  point  is  calculated.  Then  the  first,  second,  third,  or 
fourth  checkpoint  bit  is  calculated  and  stored  in  local  variable 
I0P.  Logical  checkpoint  flag  ZCKPTl  is  then  set  TRUE  if  IDP 
is  equal  to  1 .  A  test  is  then  performed  to  determine  if  debug 
output  is  desired  and  control  returns  to  the  calling  program. 
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P.5  ERROR  MESSAGE  LIST 

The  following  is  a  error  message  list  for  PMAP.  To  facilitate  error 
lookup,  the  errors  are  categorized  by  the  subroutine  name  and  the  state¬ 
ment  number  in  the  subroutine  where  the  error  was  detected.  Error  detec¬ 
tion  sometimes  generates  several  messages  providing  a  trace  to  the  source 
of  the  error.  Some  of  the  statement  numbers  were  not  properly  initialized 
prior  to  error  message  printing.  These  statement  numbers  are  listed  or 
'XXXX'  . 
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P.6  OUTPUT  RECORD  FORMATS 

PMAP  produces  a  printed  output  report  and  map  plot  and  if 
requested,  a  tape  consisting  of  the  sites  which  engage  at  least  one  path. 

P.6.1  Output  Report  and  Map  Plot 

PMAP  produces  a  site  data  table,  site  plot,  path  data  table, 
path  point  plot,  and  engageabi 1 i ty  table.  Annotated  examples  of  each  of 
these  types  of  output  may  be  found  in  Volume  111  D.  TERAIN,  PMAP,  S0RTEV 
User/Planner  Manual,  Paragraph  P.i;. 

P.6.2  Site  Tape 

If  requested  by  the  analyst,  PMAP  will  punch  one  seven 
word  logical  record  for  each  site  which  has  at  least  one  engageable  path. 

As  shown  in  Figure  P.6-1,  this  record  contains  the  system  identifier, 
site  identifier,  site  location  in  UTM  coordinates,  site  altitude,  and  the 
si te  sector. 
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P.7  MACHINE  ENVIRONMENT 

P.7*l  Hardware  Requirements 

The  MICOM  version  of  TACOS  II  will  run  on  any  computer  in  the 
CDC  6000,  7000,  or  CYBER-70  series. 

P.7.1.1  Core  Memory 

PMAP  requires  56175  (octal)  words  of  core  memory  for  execu¬ 
tion.  The  size  of  the  loader  at  the  installation  under  consideration  must 
be  accounted  for  to  obtain  the  core  requirement  for  loading.  ECS  is  not 
requi red. 

P.7.1.2  Disk  Tape  Space 

A  minimum  of  two  7“track  tape  units  are  required  during  the 
execution  of  PMAP.  The  precise  number  of  PRU's  of  disk  storage  use  will 
depend  on  the  size  of  the  terrain  file  used.  Each  100  KM  grid  square  con¬ 
tains  22  1717  word  records  and  2  1515  word  records  or  a  total  of  40,804  words 
per  grid  square.  Grid  zone  32UF  contains  12  100  KM  grid  squares,  489,648 
words  or  7,651  PRUs. 

P.7.1.3  Other  1/0  Devices 

A  card  reader  and  line  printer  are  necessary  for  normal 
1/0  interfacing  with  the  user. 

P.7-2  Operating  System  Requirements 

P.7.2.1  Basic  System 

The  operating  system  for  the  M1C0M  version  of  TACOS  II  is 
SCOPE  3.4.  The  source  language  for  PMAP  is  compiled  under  FORTRAN 
EXTENDED  V4. 

P.7.2.2  Non-Standard  Options 

PMAP  requires  no  non-standa-d  options. 
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P.7.2.3  Operator  Instructions 

If  a  non-standard  line  spacing  Is  specified,  that  is,  other 
than  6  lines  per  inch,  the  operator  should  be  informed  prior  to  program 
execution. 


1 


P.7-2 


x85< 


BRADDOCK.  DUNN  AMD  McDONALD. INC. 


APPENDIX  A 
GLOSSARY  FOR  PMAP 
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LNBUF(lA)  ZZCORD  INTEGER  Line  buffer  used  for  storing 
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